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2.0 INTRODUCTION 


A two element controller was designed for each axis to satisfy the two distinct sets of 
requirements. For the feedback loop, an integral LQG methodology [3] was used in developing 
the design to comply with feedback requirements. For the feedforward loop, an ideal model 
was developed using first or second order filters to achieve the necessary frequency separation 
between feedback and feedforward loops, and provide the desired command limiting as 
dictated by passenger comfort and airplane dynamics constraints. 

With regard to the longitudinal axis, some important differences between the design presented 
here (i.e.. Total Energy Control System - TECS - [4, 5, 6]) and the conventional control systems 
designs currently used on Boeing commercial airplanes are summarized in Table 1. The design 
presented herein is identified as BCA RESEARCH in the Table. Item definitions in the Table 
that are not self-explanatory are defined in the following sections. (NOTE: Table 1 should be 
referred to when progressing through this document.) 

With regard to the lateral axis, the design process presented here represents a departure from 
most previous techniques in that the roll and yaw axis controllers are designed simultaneously. 
Most previous techniques have traditionally developed the yaw damper first to control the 
rudder, and then designed the autopilot to control the ailerons. 


Table 1. Comparison of Controller Design Techniques 


ITEM 

BCA RESEARCH 

TECS 

CONVENTIONAL 

INTEGRATED 
PITCH / THRUST 

YES 

YES 

NO 

MODE SWITCHING 

FEEDFORWARD LOOP 

FEEDBACK LOOP 

FEEDBACK/FEEDFORWARD LOOP 

LIMITING 

FEEDFORWARD LOOP 

FEEDBACK LOOP 

FEEDBACK/FEEDFORWARD LOOP 

INITIALIZATION 

ALL STATES TO 0 

CURRENT STATE OF 
AIRPLANE 

CURRENT STATE OF 
AIRPLANE 

DYNAMIC RESPONSE 

MODE INDEPENDENT 

MODE DEPENDENT 

MODE DEPENDENT 

INTEGRAL GAINS 
(outer loop) 

HIGH 

LOW 

LOW 

PROPORTIONAL GAINS 
(inner loop) 

LOW 

HIGH 

HIGH 


PAf,f / 


y HLAN* 


2 





































1.0 SUMMARY 


The objective of this task is to design integrated autopilot controllers for both the longitudinal 
and lateral axes of the NASA TSRV airplane such that the longitudinal axis for 
autopilot/autothrottle provides independent control of altitude and speed, and the lateral axis 
for autopilot/yaw damper/sideslip controller provides independent control of heading and 
sideslip while augmenting dutch roll stability. 

The approach used to achieve this objective is the two-degree-of-freedom design philosophy 
first advocated by I. Horowitz [1]. With this philosophy, one first designs the feedback 
controller to satisfy stability and performance robustness requirements (i.e., throughout the 
flight envelope), after which one then designs the feedforward controller to satisfy 
performance requirements (Milspec [2], or other ideal airplane response specifications). Figure 
1 illustrates this concept with the arrows indicating the two degrees of freedom that the 
designer must follow to obtain a satisfactory design. 



Figure 1. Two - Degree - of - Freedom Design 


t. 





3.0 SYMBOLS AND ABBREVIATIONS 

3.1 SYMBOLS 


a 

% 

8 e 

8th 

0 

x 

y 

Yc 

l 

AFTLIM 

ALTH 

DFM 

g 

FAFLR 

FAFPA 

FAGA 

FAGLDA 

FWDLIM 

H 

H 

H c 

H c 

Hc g 

Heg 

H e 

HGCG 

HLIM 

hlim 

hrad 

11 

12 
IH 
Iv 


- Angle of attack 

- Damping ratio 

- Incremental elevator position 

- Incremental throttle position 

- Pitch altitude 

- Time constant 

- Flight path angle 

- Flight path angle command 

- Moment arm 

- Throttle aft limit flag 

- Altitude hold mode flag 

- Flap position of airplane 

- Gravity constant 

- Flare mode flag 

- FPA mode flag 

- Go-around mode flag 

- Glideslope mode arm flag 

- Throttle forward limit flag 

- Altitude 

- Altitude rate 

- Altitude command 

- Altitude rate command 

- Altitude at center of gravity for airplane 

- Altitude rate at center of gravity for airplane 

- Altitude error 

- Altitude measured at center of gravity for airplane 

- Altitude rate limiter 

- Altitude acceleration limiter 

- Radar altitude 

- Integrator #1 

- Integrator #2 

- Altitude integrator 

- Airpseed integrator 
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KTHR - Thrust to throttle handle gain ratio 
Nzcg - N z sensed center of gravity for airplane 

Q - Pitch rate 

Q - Dynamic pressure 

QC - Dynamic pressure 

s - Laplace transform 

SPDH - Speed hold mode flag 

AT - Sampling period 

u * Inertial speed 

V - True airspeed 

V - True acceleration 

V c - Airspeed command 

V c - Acceleration command 

VF - Filtered airspeed 

Vi - Inertial acceleration 

Vy jm - Acceleration command limit 

VP - True airspeed 

W - Weight of airplane 

W n - Natural frequency 
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3.2 ABBREVIATIONS 


BCA 

CAS 

eg- 

EPR 

4D 

FPA 

GSE 

IRU 

LQG 

LQR 

MCP 

MPAC 

rhp 

SOE 

TAS 

TSRV 

VTAS 


- Boeing Commercial Airplanes 

- Calibrated airspeed 

- Center of gravity 

- Engine pressure ratio 

- Four-dimensional 

- Flight path angle 

- Glideslope error 

- Inertial reference unit 

- Linear quadratic gaussian 

- Linear quadratic regulator 

- Mode control panel 

- Control analysis and design package developed at Boeing 

- Right half plane 

- Speed on elevator 

- True airspeed 

- Transport systems research vehicle 

- True airspeed 
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4.0 SETTING UP FLIGHT CONDITIONS AND AERO MODELS 

This section contains flight conditions and the procedure by which these conditions were 
developed for this study. The first step in setting up flight conditions was to create a fairly 
accurate definition of the flight envelope from the operational manual of the 737-200 airplane. 
Tables 2 and 3 were generated from the data shown in Figure 2 (stall speeds) and Figure 3 
(operational limits and placards). 

Table 2. Light Weight - 80,000 lbs. with C. G. Variation of .05 to .3 


FLAPS 

SETTING 

V stall 

PLACARD 

SPEED 

GEAR 

0 

128 

— 

UP 

1 

106 

210 

UP 

5 

100 

210 

UP 

10 

96 

210 

UP 

15 

94 

195 

UP 

25 

91 

190 

UP 

30 

89 

185 

DN 

40 

87 

170 

DN 


Table 3. Heavy Weight - 110,000 lbs. with C. G. Variation of .05 to .3 


FLAPS 

SETTING 

v stall 

PLACARD 

SPEED 

GEAR 

0 

150 

— 

UP 

1 

125 

210 

UP 

5 

118 

210 

UP 

10 

114 

210 

UP 

15 

110 

195 

UP 

25 

108 

190 

UP 

30 

106 

185 

DN 

40 

103 

170 

DN 
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Figure 3. Operational Limits and Placards 










Using Tables 2 and 3, the operational envelope was defined by the flight condition files for: 


a . Cruise 

b. Cruise Transition 

c. Landing 

4.1 CRUISE 

The cruise flight conditions file was constructed using the following data (from Tables 2 and 3: 


Flaps 

0 

Altitudes - 

10000 ft, 25000 ft and 35000 ft 

Weights 

maximum take-off weight and minimum in-flight weight plus 10000 lbs 

C. G.'s 

maximum forward and aft C. G.s (i.e., .05 and .3 respectively) 

Gear 

landing gear up 

Speeds 

The low speed limit for 1.3 g maneuver margin to initial buffet, maximum 
operating speed (M mo or V mo ) or highest attainable Mach number, and 
two intermediate speeds (chosen such that there is equal Mach number 
spacing between the four Mach numbers). 

4.2 CRUISE TRANSITION 

The cruise transition flight conditions file was constructed using the following data along with 
Tables 2 and 3: 

Flaps 

1, 5, 10, 15, 25 

Altitude 

0 ft barometric altitude 

Weights 

maximum landing weight plus 10000 lbs and minimum in-flight weight 
plus 10000 lbs 

C. G.s 

maximum forward and aft C. G.s (i.e., .05 and .3 respectively) 

Gear 

landing gear up 

Speeds 

1.3 V s (V s = stall speed), Vf p = Flap Placard Speed, and two 
intermediate speeds (1.3 V s and Vfp-chosen such that there is 
approximately equal Qc spacing between the four speeds). 
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4.3 LANDING 


The landing flight conditions of Glideslope and Flare were constructed using the following data 
along with Tables 2 and 3: 

Flaps - 30, 40 

Altitude - 0 ft. barometric altitude 

Weights - maximum landing weight plus 10000 lbs and minimum in-flight weight 
plus 10000 lbs 

C.G.s - maximum forward and aft C. G.s 

Gear - landing gear down 

Speeds - V app> v app + 5 kts, V app + 20 kts (V app = approach speed) 

The definition and limitation of the parameters for these flight conditions varied slightly due 
to the trimming of the airplane on the Harris Simulation Computer. Table 4a-c show the 
detailed definitions of 48 flight conditions for cruise, 80 flight conditions for cruise transition, 
and 32 flight conditions for landing. 

4.4 GENERATION OF THE OPEN-LOOP AIRPLANE MODELS 


Each of the flight conditions was read to the Harris Simulator Computer and the program 
MATGEN (matrix generation) executed to obtain the A and B matrices of the linear airplane 
models. The A and B matrices were formed by perturbing predefined states and controls when 
MATGEN was run on the Harris Simulator. 

The procedure for MATGEN on the Harris is: 

a. Aircraft is trimmed at specified flight condition. 

b. Perturb state or surface a small amount. 

c. Allow accelerations to settle. 

d. Store resultant changes in the airplane states. 

e. Develop an A matrix for each perturbation (i.e., positive perturbation, negative 
perturbation, and average of the two). 

f . Develop B matrix similarly to A for each control surface or discrete perturbed. 
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Table 4a. Cruise Flight Conditions 


COND. 

MACH 

VEq 

Qe 

VTA S 

ALT. 

CG. 

WEIGHT 

GEAR 

FLAPS 

COND56789012345678 9012345678 901234567890123456789012345678 9012345678 901234567 8 90123456789012345678 90 

1 

.30999 

170.67 

100.27 

333.98 

10000. 

5.0000E-02 

80000. 

0 . 

0 . 

2 

.41999 

231.92 

187.76 

452.50 

10000. 

5.0000E-02 

80000. 

0 . 

0 . 

3 

.52999 

293.75 

306.82 

571.01 

10000. 

5.0000E-02 

80000. 

0 . 

0 . 

4 

.63999 

356.26 

461.76 

689.52 

10000. 

5.0000E-02 

80000. 

0 . 

0 . 

5 

.30999 

170.67 

100.27 

333.98 

10000. 

.3 

80000. 

0 . 

0 . 

6 

.41999 

231.92 

187.76 

452.50 

10000. 

.3 

80000. 

0 . 

0 . 

7 

.52999 

293.75 

306.82 

571.01 

10000. 

.3 

80000. 

0 . 

0 . 

8 

.63999 

356.26 

461.76 

689.52 

10000. 

.3 

80000. 

0 . 

0 . 

9 

.36999 

204.01 

144.30 

398.63 

10000. 

5.0000E-02 

110000. 

0 . 

0 . 

10 

.45999 

254.33 

227.20 

495.59 

10000. 

5.0000E-02 

110000. 

0 . 

0 . 

11 

.54999 

305.06 

332.18 

592.56 

10000. 

5 . 0000E-O2 

110000. 

0 . 

0 . 

12 

.63999 

356.26 

461.76 

689.52 

10000. 

5.0000E-02 

110000. 

0 . 

0 . 

13 

.36999 

204.01 

144.30 

398.63 

10000. 

.3 

110000. 

0 . 

0 . 

14 

.45999 

254.33 

227.20 

495.59 

10000. 

.3 

110000. 

0 . 

0 . 

15 

.54999 

305.06 

332.18 

592.56 

10000. 

.3 

110000. 

0 . 

0 . 

16 

.63999 

356.26 

461.76 

689.52 

10000. 

.3 

110000. 

0 . 

0 . 

17 

.47999 

196.92 

134.22 

487.66 

25000. 

5.0000E-02 

80000. 

0 . 

0 . 

18 

.53999 

222.53 

172.45 

548.62 

25000. 

5.00.00E— 02 

80000. 

0 . 

0 . 

19 

.59999 

248.48 

216.52 

609.58 

25000. 

5.0000E-02 

80000. 

0 . 

0 . 

20 

.83999 

356.18 

461.56 

853.41 

25000. 

5.0000E-02 

80000. 

0 . 

0 . 

21 

.47999 

196.92 

134.22 

487.66 

25000. 

.3 

80000. 

0 . 

0 . 

22 

.53999 

222.53 

172.45 

548.62 

25000. 

.3 

80000. 

0 . 

0 . 

23 

.59999 

248.48 

216.52 

609.58 

25000. 

.3 

80000. 

0 . 

0 . 

24 

.83999 

356.18 

461.56 

853.41 

25000. 

.3 

80000. 

0 . 

0 . 

25 

.55999 

231.14 

186.47 

568.94 

25000. 

5.0000E-02 

110000. 

0 . 

0 . 

26 

.65999 

274.80 

266.88 

670.54 

25000. 

5.0000E-02 

110000. 

0 . 

0 . 

27 

.74999 

315.01 

355.45 

761.98 

25000. 

5.0000E-02 

110000. 

0 . 

0 . 

28 

.83999 

356.18 

461.56 

853.41 

25000. 

5.0000E-02 

110000. 

0 . 

0 . 

29 

.55999 

231.14 

186.47 

568.94 

25000. 

.3 

110000. 

0 . 

0 . 

30 

.65999 

274.80 

266.88 

670.54 

25000. 

.3 

110000. 

0 . 

0 . 

31 

.74999 

315.01 

355.45 

761.98 

25000. 

.3 

110000. 

0 . 

0 . 

32 

.83999 

356.18 

461.56 

853.41 

25000. 

.3 

110000. 

0 . 

0 . 

33 

.54999 

181.59 

113.76 

535.08 

35000. 

5.0000E-02 

80000. 

0 . 

0 . 

34 

.64999 

216.92 

163.63 

632.37 

35000. 

5 . OOOOE-02 

80000. 

0 . 

0 . 

35 

.74999 

253.36 

225.42 

729.66 

35000. 

5.0000E-02 

80000. 

0 . 

0 . 

36 

.83999 

287.22 

292.71 

817.22 

35000. 

5.0000E-02 

80000. 

0 . 

0 . 

37 

.54999 

181.59 

113.76 

535.08 

35000. 

.3 

80000. 

0 . 

0 . 

38 

.64999 

216.92 

163.63 

632.37 

35000. 

.3 

80000. 

0 . 

0 . 

39 

.74999 

253.36 

225.42 

729.66 

35000. 

.3 

80000. 

0 . 

0 . 

40 

.83999 

287.22 

292.71 

817.22 

35000. 

.3 

80000. 

0 . 

0 . 

41 

.65999 

220.51 

169.25 

642.10 

35000. 

5.0000E-02 

110000. 

0 . 

0 . 

42 

.71999 

242.30 

205.53 

700.47 

35000. 

5.0000E-02 

110000. 

0 . 

0 . 

43 

.77999 

264.53 

246.54 

758.85 

35000. 

5.0000E-02 

110000. 

0 . 

0 . 

44 

.83999 

287.22 

292.71 

817.22 

35000. 

5.0000E-02 

110000. 

0 . 

0 . 

45 

.65999 

220.51 

169.25 

642.10 

35000. 

.3 

110000. 

0 . 

0 . 

46 

.71999 

242.30 

205.53 

700.47 

35000. 

.3 

110000. 

0 . 

0 . 

47 

.77999 

264.53 

246.54 

758.85 

35000. 

.3 

110000. 

0 . 

0 . 

48 

.83999 

287.22 

292.71 

817.22 

35000. 

.3 

110000. 

0 . 

0 . 


Table 4b. Cruise Transition Flight Conditions 


COND. 

MACH 

VEq 

Qe 

VTAS 

ALT. 

CG. 

WEIGHT 

GEAR FLAPS 

COND567 8 9012345 678 9012345 678 9012345678 901234567 89012345678 90123456789012345678 90123456789012345678 90 


1 

.20862 

138.00 

65.180 

232.92 

100. 

5. OOOOE—02 

80000. 

1 . 

0 . 

2 

.25095 

166.00 

94.771 

280.17 

100. 

5 . 000 0E— 02 

80000. 

1 . 

0 . 

3 

.28572 

189.00 

123.42 

318.99 

100. 

5. 0000E-02 

80000. 

1 . 

0 . 

4 

.31747 

210.00 

153.10 

354.44 

100. 

5. OOOOE-02 

80000. 

1 . 

0 . 

5 

.20862 

138.00 

65.180 

232.92 

100. 

.3 

80000. 

1 . 

0 . 

6 

.25095 

166.00 

94.771 

280.17 

100. 

.3 

80000. 

1 . 

0 . 

7 

.28572 

189.00 

123.42 

318.99 

100. 

.3 

80000. 

1 . 

0 . 

8 

.31747 

210.00 

153.10 

354.44 

100. 

.3 

80000. 

1 . 

0 . 

9 

.24642 

163.00 

91.325 

275.11 

100. 

5. OOOOE-02 

110000. 

1 . 

0 . 

10 

.27212 

180.00 

111.73 

303.80 

100. 

5. OOOOE-02 

110000. 

1 . 

0 . 

11 

.29328 

194.00 

130.18 

327.43 

100. 

5. OOOOE-02 

110000. 

1 . 

0 . 

12 

.31747 

210.00 

153.10 

354.44 

100. 

5.0000E— 02 

110000. 

1 . 

0 . 

13 

.24642 

163.00 

91.325 

275.11 

100. 

.3 

110000. 

1 . 

0 . 

14 

.27212 

180.00 

111.73 

303.80 

100. 

.3 

110000. 

1 . 

0 . 

15 

.29328 

194.00 

130.18 

327.43 

100. 

.3 

110000. 

1 . 

0 . 

16 

.31747 

210.00 

153.10 

354.44 

100. 

.3 

110000. 

1 . 

0 . 

17 

.19653 

130.00 

57.771 

219.41 

100. 

5. OOOOE-02 

80000. 

5. 

0 . 

18 

.24339 

161.00 

89.065 

271.74 

100. 

5. OOOOE-02 

80000. 

5. 

0 . 

19 

.28421 

188.00 

122.09 

317.31 

100. 

5. 0000E-02 

80000. 

5. 

0 . 

20 

.31747 

210.00 

153.10 

354.44 

100. 

5. OOOOE-02 

80000. 

5. 

0 . 

21 

.19653 

130.00 

57.771 

219.41 

100. 

.3 

80000. 

5. 

0 . 

22 

.24339 

161.00 

89.065 

271.74 

100. 

.3 

80000. 

5. 

0 . 

23 

.28421 

188.00 

122.09 

317.31 

100. 

.3 

80000. 

5. 

0 . 

24 

.31747 

210.00 

153.10 

354.44 

100. 

.3 

80000. 

5. 

0 . 

25 

.23130 

153.00 

80.319 

258.23 

100. 

5 . 0000E-02 

110000. 

5. 

0 . 

26 

.26304 

174.00 

104.28 

293.68 

100. 

5. OOOOE-02 

110000. 

5. 

0 . 

27 

.29328 

194.00 

130.18 

327.43 

100. 

5. 0000E-02 

110000. 

5. 

0 . 

28 

.31747 

210.00 

153.10 

354.44 

100. 

5. OOOOE-02 

110000. 

5. 

0 . 

29 

.23130 

153.00 

80.319 

258.23 

100. 

.3 

110000. 

5. 

0 . 

30 

.26304 

174.00 

104.28 

293.68 

100. 

.3 

110000. 

5. 

0 . 

31 

.29328 

194.00 

130.18 

327.43 

100. 

.3 

110000. 

5. 

0 . 

32 

.31747 

210.00 

153.10 

354.44 

100. 

.3 

110000. 

5. 

0 . 

33 

.18897 

125.00 

53.374 

210.97 

100. 

5. 0000E-02 

80000. 

10. 

0 . 

34 

.24037 

159.00 

86.834 

268.36 

100. 

5. OOOOE-02 

80000. 

10. 

0 . 

35 

.28270 

187.00 

120.77 

315.62 

100. 

5. 0000E-02 

80000. 

10. 

0 . 

36 

.31747 

210.00 

153.10 

354.44 

100. 

5. OOOOE-02 

80000. 

10. 

0 . 

37 

.18897 

125.00 

53.374 

210.97 

100. 

.3 

80000. 

10. 

0 . 

38 

.24037 

159.00 

86.834 

268.36 

100. 

.3 

80000. 

10. 

0 . 

39 

.28270 

187.00 

120.77 

315.62 

100. 

.3 

80000. 

10. 

0 . 

40 

.31747 

210.00 

153.10 

354.44 

100. 

.3 

80000. 

10. 

0 . 

41 

.22374 

148.00 

75.091 

249.79 

100. 

5. OOOOE-02 

110000. 

10. 

0 . 

42 

.26002 

172.00 

101.86 

290.30 

100. 

5. 0000E-02 

110000. 

10. 

0 . 

43 

.29026 

192.00 

127.45 

324.06 

100. 

5. 0000E-02 

110000. 

10. 

0 . 

44 

.31747 

210.00 

153.10 

354.44 

100. 

5. 0000E-02 

110000. 

10. 

0 . 

45 

.22374 

148.00 

75.091 

249.79 

100. 

.3 

110000. 

10. 

0 . 

46 

.26002 

172.00 

101.86 

290.30 

100. 

.3 

110000. 

10. 

0 . 

47 

.29026 

192.00 

127.45 

324.06 

100. 

.3 

110000. 

10. 

0 . 

48 

.31747 

210.00 

153.10 

354.44 

100. 

.3 

110000. 

10. 

0 . 

49 

.18443 

122.00 

50.821 

205.91 

100. 

5. 0000E-02 

80000. 

15. 

0 . 

50 

.22979 

152.00 

79.259 

256.54 

100. 

5. OOOOE-02 

80000. 

15. 

0 . 

51 

.26607 

176.00 

106.74 

297.05 

100. 

5. OOOOE-02 

80000. 

15. 

0 . 

52 

.29479 

195.00 

131.55 

329.12 

100. 

5. 0000E-02 

80000. 

15. 

0 . 

53 

.18443 

122.00 

50.821 

205.91 

100. 

.3 

80000. 

15. 

0 . 

54 

.22979 

152.00 

79.259 

256.54 

100. 

.3 

80000. 

15. 

0 . 

55 

.26607 

176.00 

106.74 

297.05 

100. 

.3 

80000. 

15. 

0 . 

56 

.29479 

195.00 

131.55 

329.12 

100. 

.3 

80000. 

15. 

0 . 

57 

.21618 

143.00 

70.044 

241.35 

100. 

5. 0000E-02 

110000. 

15. 

0 . 

58 

.24339 

161.00 

89.065 

271.74 

100. 

5. 0000E-02 

110000. 

15. 

0 . 

59 

.27212 

180.00 

111.73 

303.80 

100. 

5. 0000E-02 

110000. 

15. 

0 . 

60 

.29479 

195.00 

131.55 

329.12 

100. 

5. 0000E-02 

110000. 

15. 

0 . 

61 

.21618 

143.00 

70.044 

241.35 

100. 

.3 

110000. 

15. 

0 . 

62 

.24339 

161.00 

89.065 

271.74 

100. 

.3 

110000. 

15. 

0 . 

63 

.27212 

180.00 

111.73 

303.80 

100. 

.3 

110000. 

15. 

0 . 

64 

.29479 

195.00 

131.55 

329.12 

100. 

.3 

110000. 

15. 

0 . 

65 

.17839 

118.00 

47.517 

199.16 

100. 

5 . 0000E-02 

80000. 

25. 

0 . 

66 

.22223 

147.00 

74.067 

248.11 

100. 

5. 0000E-02 

80000. 

25. 

0 . 

67 

.25700 

170.00 

99.470 

286.93 

100. 

5. OOOOE-02 

80000. 

25. 

0 . 

68 

.28723 

190.00 

124.76 

320.68 

100. 

5. OOOOE-02 

80000. 

25. 

0 . 

69 

.17839 

118.00 

47.517 

199.16 

100. 

.3 

80000. 

25. 

0 . 

70 

.22223 

147.00 

74.067 

248.11 

100. 

.3 

80000. 

25. 

0 . 

71 

.25700 

170.00 

99.470 

286.93 

100. 

.3 

80000. 

25. 

0 . 

72 

.28723 

190.00 

124.76 

320.68 

100. 

.3 

80000. 

25. 

0 . 

73 

.21164 

140.00 

67.104 

236.29 

100. 

5. OOOOE-02 

110000. 

25. 

0 . 

74 

.24037 

159.00 

86.834 

268.36 

100. 

5. OOOOE-02 

110000. 

25. 

0 . 

75 

.26607 

176.00 

106.74 

297.05 

100. 

5 . OOOOE-02 

110000. 

25. 

0 . 

76 

.28723 

190.00 

124.76 

320.66 

100. 

5. OOOOE-02 

110000. 

25. 

0 . 

77 

.21164 

140.00 

67.104 

236.29 

100. 

.3 

110000. 

25. 

0 . 

78 

.24037 

159.00 

86.834 

268.36 

100. 

.3 

110000. 

25. 

0 . 

79 

.26607 

176.00 

106.74 

297.05 

100. 

.3 

110000. 

25. 

0 . 

80 

.28723 

190.00 

124.76 

320.68 

100. 

.3 

110000. 

25. 

0 . 
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Table 4c. Landing Flight Conditions 


COND. MACH 

VEq 

Qb 

VTAS 

ALT. 

CG. 

WEIGHT 

GEAR FLAPS 

COND5 67 8 90 12345 678 9012345 678 901234567 8 9012345678 9012345678 9012345678 9012345 67 8 90123 4 567 8 90 12 3 45 678 90 

1 

.17536 

116.00 

45.908 

195.78 

100. 

5.Q000E-02 

80000. 

30. 

1. 

2 

.21618 

143.00 

70.044 

241.35 

100. 

5. 0000E— 02 

80000. 

30. 

1. 

3 

.24793 

164.00 

92.467 

276.80 

100. 

5.0000E-02 

80000. 

30. 

1. 

4 

.27967 

185.00 

118.15 

312.24 

100. 

5.0000E-02 

80000. 

30. 

1. 

5 

.17536 

116.00 

45.908 

195.78 

100. 

.3 

80000. 

30. 

1. 

6 

.21618 

143.00 

70.044 

241.35 

100. 

.3 

80000. 

30. 

1. 

7 

.24793 

164.00 

92.467 

276.80 

100. 

.3 

80000. 

30. 

1. 

8 

.27967 

185.00 

118.15 

312.24 

100. 

-3 

80000. 

30. 

1. 

9 

.20862 

138.00 

65.180 

232.92 

100. 

5.0000E-02 

110000. 

30. 

1. 

10 

.23432 

155.00 

82.461 

261.61 

100. 

5 . OOOOE— 02 

110000. 

30. 

1. 

11 

.25700 

170.00 

99.470 

286.93 

100. 

5.0000E-02 

110000. 

30. 

1. 

12 

.27967 

185.00 

118.15 

312.24 

100. 

5 . OOOOE— 02 

110000. 

30. 

1. 

13 

.20862 

138.00 

65.180 

232.92 

100. 

.3 

110000. 

30. 

1. 

14 

.23432 

155.00 

82.461 

261.61 

100. 

.3 

110000. 

30. 

1. 

15 

.25700 

170.00 

99.470 

286.93 

100. 

.3 

110000. 

30. 

1. 

16 

,27967 

185.00 

118.15 

312.24 

100. 

.3 

110000. 

30. 

1. 

17 

.17083 

113.00 

43.547 

190.72 

100. 

5 . OOOOE— 02 

80000. 

39.99 

1. 

18 

.20257 

134.00 

61.418 

226.16 

100. 

5 . OOOOE— 02 

80000. 

39.99 

1. 

19 

.23281 

154.00 

81.386 

259.92 

100. 

5 . OOOOE— 02 

80000. 

39.99 

1. 

20 

.25700 

170. Ou 

99.470 

286.93 

100. 

5 . OOOOE— 02 

80000. 

39.99 

1. 

21 

.17083 

113.00 

43.547 

190.72 

100. 

.3 

80000. 

39.99 

1. 

22 

.20257 

134.00 

61.418 

226.16 

100. 

.3 

80000. 

39.99 

1. 

23 

.23281 

154.00 

81.386 

259.92 

100. 

.3 

80000. 

39.99 

1. 

24 

.25700 

170.00 

99.470 

286.93 

100. 

.3 

80000. 

39.99 

1. 

25 

.20257 

134.00 

61.418 

226.16 

100. 

5 .OOOOE— 02 

110000. 

39.99 

1. 

26 

.22525 

149.00 

76.122 

251.48 

100. 

5.0000E-02 

110000. 

39.99 

1. 

27 

.24037 

159.00 

86.834 

268.36 

100. 

5.0000E-02 

110000. 

39.99 

1. 

28 

.25700 

170.00 

99.470 

286.93 

100. 

5.0000E-02 

110000. 

39.99 

1. 

29 

.20257 

134.00 

61.418 

226.16 

100. 

.3 

110000. 

39.99 

1. 

30 

.22525 

149.00 

76.122 

251.48 

100. 

.3 

110000. 

39.99 

1. 

31 

.24037 

159.00 

86.834 

268.36 

100. 

.3 

110000. 

39.99 

1. 

32 

.25700 

170.00 

99.470 

286.93 

100. 

.3 

110000. 

39.99 

1^ 
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All states are perturbed +1 unit, except UA which is perturbed +5 ft/sec. All surfaces are 
perturbed +1 unit, except speed brake handle which is moved +.5 and engine thrust which is 
changed via throttle handle with movements of ±100 lb. 

The resultant A and B matrices are: 


A MATRIX 


UaQe P P ♦ R 



B MATRIX 


DE DSPL DTH AIL RUD 



The A and B matrices are then expanded to include additional states and wind disturbances and 
the output matrices, C and D, are formed. 

For the longitudinal axis, the altitude state,H, is added to the A matrix and gust disturbances, 
Ug and otg,are added to the B matrix. For the lateral directional axis the yaw and lateral 
position deviation states, \j / and y, are formed and added to the A matrix, and Pg is added to 
the B matrix; the outputs p, Ny, R and yGT are formed and used to derive the C and D matrices. 
The final ABCD matrices for each axis of each flight condition are then printed. Figures 2 and 
3 show a sample case of longitudinal axis ABCD matrices and lateral directional ABCD 
matrices respectively. 
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5.0 CONTROL DESIGN MEraODOLOGY 


This section presents the control system design methodology used for this task. This design 
methodology has been developed jointly by the Boeing Advanced Systems (BAS) flight controls 
research and the Boeing Commercial Airplanes (BCA) guidance and control research groups. 
Sections 7 and 8 document the application of this design approach to the longitudinal and 
lateral axes of the NASA TRSV airplane respectively. A model emulating the integral linear 
quadratic guassian control design technique is used in conjunction with a systematic top down 
approach. Strong emphasis is placed on understanding the dynamics of the open loop airplane 
(including coupling, controllability, and observability) and the control task requirements in 
order to develop specific design requirements that are responsive to the guidance and control 
problems without overly constraining the problem or violating the controllability limits of the 
open loop system. The following paragraphs highlight the key ingredients of this Boeing 
control system design methodology which are presented in the order that these ingredients are 
used during the design process. 

5.1 SPECIFY DESIGN OBJECTIVES 

The first step in designing of a control system is to specify the top level design objectives. These 
start with the vehicle mission goals. For example, the number of passengers, the range of 
operation, the range of the flight envelope, and the relative importance of speed, fuel economy, 
ride comfort, and safety. 

Although the controls engineer may not have any input in the initial development of the 
mission goals, it is important that these objectives be understood. All successive requirements, 
both those imposed upon and those developed by the controls engineer, must be responsive to the 
top level mission goals. Specifically, the definition of the airplane configuration and the 
architecture of the control system must be in keeping with the vehicle mission goals. In some 
instances, it may be the responsibility of the controls engineer to indicate whether the vehicle 
configuration or the the control system architecture may impede in achieving the top level 
mission goals that the control system design cannot overcome alone. When this is the case, a 
thorough understanding of the airplane and its mission goals will allow the engineer to make 
such a determination. 
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The goals for the control system are derived directly from the vehicle mission goals. An 
example of this is automatic landing in which modem commercial aircraft must operate in 
conditions where normal piloted landings are not possible because the pilot cannot see the 
ground to judge the approach. In order to continue operation in severely reduced visibility 
conditions, the plane must be able to land on its own. Therefore the requirements for an 
automatic landing system follow directly for the mission goal to be able to operate in dense fog. 

The pilot/ airplane interface must be defined once the top level control system goals have been 
determined since the mission and control system goals place demands on the airplane control 
system. The pilot, by the nature of his job, is an integral part of the system that controls the 
airplane. The control task must be divided between the pilot and the flight control system. 

The days when there was a sharp delineation between manual flight and autopilot flight on 
commercial transport aircraft are gone. Today the control systems on modern commercial 
aircraft have active stability augmentation systems, in the form of yaw dampers, operating 
during the entire flight (including manually piloted flight). Up front definition of the 
pilot/airplane interface is essential since the controller design process steps, as described in the 
following paragraphs, are highly dependent on the pilot/ airplane interface definition. 

A significant feature of the design methodology is separation of the feedforward and feedback 
controllers. Integral feedback is used to decouple steady state responses and provide the needed 
bandwidth and damping. The feedforward controller is used to decouple and shape the 
transient responses. The design task is split into two subtasks. First, the feedback controller is 
designed to achieve the desired stability characteristics, and then the feedforward controller 
is designed. The stability characteristics of the closed loop system are not affected by 
feedforward controller modifications because of the separation inherent in the controller 
structure. Within the bandwidth limitations of the feedback controller, the transient response 
and performance characteristics of the augmented airplane can be modified without affecting 
the stability and disturbance rejection characteristics. 

52 FORMULATE DESIGN REQUIREMENTS 

Once the design objectives are defined, the next step is formulation of the specific design 
requirements. The requirements must be responsive to the design objectives, as described in 
Section 2, without being overly constraining. Each design requirement must be traceable back to 
the design objective that motivated it. 
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Separation of the feedforward and feedback controllers allows for a similar separation of the 
feedforward and feedback design requirements. The feedforward controller must be responsive 
to the command response characteristics while the feedback controller must meet the 
disturbance rejection and robustness requirements. 

5.3 OPEN LOOP ANALYSIS 

A complete understanding of the open loop airplane is necessary prior to designing the 
controller. Open loop analysis gives insight into the difficulty of achieving the design 
requirements before launching into the detailed design cycle. Characteristics of the open loop 
system that will render the design requirements unachievable (e.g., insufficient control 
authority) can be identified rapidly. 

The specific open loop analyses performed are: eigenstructure decomposition, controllability 
analysis, observability analysis, computation of open loop frequency responses, and open loop 
time domain simulation. Each of these analysis tools provides a different look at the 
characteristics of the open loop system. In many cases the information provided by one type of 
analysis is the same as that provided by another. The objective of open loop analysis is to 
provide sufficient perspectives of the open loop system to allow the designer on understanding 
of the system characteristics in order to design a controller. The more that is understood of the 
system prior to designing, then the more efficiently the designer will be during the design 
phase. Each of these analysis tools is described in paragraphs 5.3.1 - 5.3.5. 

5.3.1 EIGENSTRUCTURE DECOMPOSITION 

Eigenstructure decomposition reveals the open loop mode shapes and their frequencies. Many of 
the design requirements can be translated into constraints on the closed loop eigenstructure. 
Knowledge of the open loop eigenstructure and the desired closed loop eigenstructure reveals 
the magnitude of the controller task in terms of how much the eigenstructure must be changed. 

5.3.2 CONTROLABILITY ANALYSIS 

Controllability analysis consists of computing and analyzing the controllability matrix. The 
controllability matrix, as used here, is the B matrix of the modalized open loop model: 

Controllability matrix = T'l B 
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where 


T = modal transformation matrix 
B = input-to-state rate matrix 

This analysis reveals the effect of each of the open loop system inputs on the each of its modes. 
Modes whose eigenstructure cannot be modified, because none of the inputs affects them, can be 
identified. Control inputs with identical effect on the eigenstructure, and thus only a single 
degree of input freedom among them, can be identified. 

5.3.3 OBSERVABILITY ANALYSIS 

Observability analysis is the output dual of controllability analysis. The observability matrix 
is the C matrix of the modalized open loop model: 

Observability matrix = CT 

where 

C = state-to-output matrix 
T - modal transformation matrix 

The observability matrix indicates which modes are measurable via which outputs. In order 
for an open loop mode to be successfully modified via feedback control, it must be both 
observable and controllable. 

5.3.4 FREQUENCY RESPONSE 

Two types of open loop frequency responses are computed: control effector responses and 
disturbance input responses. The control effector responses show the effect of the control inputs 
on the airplane states and outputs. These analyses reveal the ability of each control input to 
control the various airplane states and outputs at different frequencies. The phase 
relationships reveal the need for lead or lag compensation when using classical design 
techniques. 

The disturbance response characteristics revealed by the second set of open loop frequency 
responses illustrate the need for active disturbance rejection. Often closed loop disturbance 
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response requirements can be met without active compensation if the open loop responses are 
sufficiently small. Active control will be required for all outputs whose open loop responses to 
disturbances fail to meet the closed loop requirements. 

5.3.5 TIME DOMAIN SIMULATION 

Open loop time domain simulation is used to understand the characteristic responses of the open 
loop system to control and disturbance inputs. Although the information revealed via 
simulation is the same as that given for frequency response analysis, some designers are more 
familiar with the time domain. 


5.4 CONTROLLER STRUCTURE 

The control structure used for this design task (see Figure 4) is a two-degree-of-freedom 
approach which maintains separation between the feedforward and feedback controllers. The 
function of the feedback controller is to provide the necessary stability augmentation, sufficient 
command response bandwidth to satisfy the performance requirements, and robustness to 
parameter uncertainties. Whereas the function of the feedforward controller is to shape the 
pilot inputs such that the closed-loop performance requirements are met. 

The feedback controller has a full state integral structure. A regulated variable is chosen for 
each independent control input (see paragraph 5.5.2). The feedback controller places an 
integrator on each of the regulated variables that is not itself the output of an integrator 
within the open loop system. The feedback gain matrix contains gains from each of the open 
loop model states and the regulated variable integrators to each of the independent control 
inputs. Feedback signal estimation is included for those signals not directly available over the 
required frequency range using sensor outputs. 

The feedforward controller consists of an ideal model defining the desired regulated variable 
response to each pilot input command. Command limiting, envelope protection, and transient 
smoothing are all part of the feedforward controller. In some cases the feedforward controller 
will send cross feeds to two different regulated variables to decouple the doseu-loop responses. 
An example of this would be to feed the heading command to sideslip and improve turn 
coordination during turn initiation. Separation of the feedforward and feedback controllers 
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High gain at low frequency will result in large penalties for steady state errors. The Riccati 
solution will provide good steady state tracking. Low gain at high frequency will result in 
small penalties for high frequency errors. The Riccati solution will ignore high frequency 
errors, thus avoiding excessive bandwidths that would lead to input actuator limiting. 

Since the synthesis model is square (i.e., the number of inputs equals the number of outputs), its 
transmission zeros can be computed. With LQG design, the transmission zeros of the synthesis 
model designate the asymptotic locations for the closed loop poles. Furthermore, zeros created 
when forming criteria outputs become transmission zeros of the synthesis model. This feature is 
used during criteria output formation to establish targets for closed loop poles whose open loop 
characteristics are not satisfactory. An example is the pair of complex zeros added to the 
sideslip criterion output to attract the dutch roll mode. 

The primary concerns during formation of the criteria outputs are the input/output frequency 
response shape (gain and phase) and the placement of transmission zeros. Selection of the 
weighting values (see paragraph 5.6) is used to scale the magnitudes of the synthesis model 
input/output frequency responses. 


5.6 COST FUNCTION WEIGHTINGS, CONTROLLER GAINS, AND GAIN SCHEDULING 

Once the synthesis model is formed, the next task is to determine the feedback gains by 
choosing the cost function weightings, solving the Riccati equation to specify the gains, and 
using gain scheduling to define the controller for a range of flight conditions. The emphasis 
with this design methodology is directed at making the tasks described herein as simple as 
possible. The bulk of the design effort is spent designing the criteria outputs. Once the 
synthesis model is complete, computation of the feedback gains is automated by the Riccati 
equation. 

5.6.1 COST FUNCTION WEIGHTINGS 

For LQG design the synthesis model is used by the Riccati equation in conjunction with diagonal 
input and output weighting matrices to minimize the following quadratic cost function: 

J = 1/2 (Q T yQ + R T uR)d 
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isolates the stability and performance characteristics of the closed loop system. Modifications 
made to the feedforward controller will not affect system stability. 

5.5 SYNTHESIS MODEL DEFINITION 


The synthesis model is used as an input to the Riccati equation which solves for the full state 
feedback gains. The synthesis model is built using the open loop model as its core. Criteria 
outputs are formed for output weighting with the Riccati equation. 

5.5.1 INDEPENDENT INPUTS 


The first step in forming the synthesis model is identification of the independent control inputs 
to be driven by the controller. It is important to verify via controllability analysis that the 
selected inputs are independent. For the longitudinal and lateral examples presented in 
Sections 7 and 8, two independent control inputs are used in each case: elevator and throttle are 
used for the longitudinal design, whereas ailerons and rudder are used for the lateral design. 

5.5.2 REGULATED VARIABLES 


Once the independent inputs are chosen, the next task is selection of regulated variables. A 
regulated variable is one that is to be controlled in steady state. One regulated variable is 
chosen for each independent input. Regulated variables must be available as sensor outputs (or 
estimates) for feedback. Of particular importance is the low frequency integrity of regulated 
variable signals since they will be integrated to provide steady state tracking. In the 
longitudinal case, the regulated variables are altitude and speed. Heading and sideslip are 
used for the lateral axis synthesis model. 

5.5.3 CRITERIA OUTPUTS 

The final step in building the synthesis model is formation of the criteria outputs. One criteria 
output is formed for each regulated variable. An integrator is added to each regulated variable 
output that is not itself the output of an integrator in the open loop system. Criteria outputs are 
then formed by adding together the regulated variable, its integrator, and any other open loop 
system outputs. Scaling between the signals that compose each criteria output are chosen with 
the objective of constructing the criteria outputs such that the frequency responses from the 
inputs to the criteria outputs exhibit high gain at low frequency and roll off at high frequency. 
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where 


y = vector of criteria outputs 
u = vector or control inputs 
Q = diagonal criteria output weighting matrix 
R = diagonal control input weighting matrix 

The diagonal Q matrix weights the criteria outputs which are composed of linear combinations 
of the synthesis model states. The same cost faction could be realized by applying the 
appropriate full Q' matrix to the synthesis model state vector. The design process used here 
simplifies the weighting matrix selection task, without loss of generality, by first forming 
criteria outputs and then using a diagonal Q rather than the full Q' along with the state vector. 

Cost function weightings are used to place the bandwidth of the closed-loop system. Increasing 
the weighting on an input will lower the bandwidth, while increasing an output weighting 
will increase the bandwidth. The degree to which changing the relative weighting magnitude 
between inputs (or outputs) changes the closed-loop system will depend on the coupling within 
the synthesis model. Relative weighting magnitude will determine the level of cross coupling 
of the closed loop response. For example, the amount of sideslip response to a heading 
command. 

Initial values for the weighting matrices are chosen by looking at the open loop frequency 
responses of the synthesis model. The closed loop bandwidths will be approximately the same 
as the bandwidths of these open loop responses. Iteration is used, looking at the closed-loop 
system characteristics, to arrive at the final set of weightings. 

5.6.2 CONTROLLER GAINS AND GAIN SCHEDULING 

Once the weighting matrices are selected the feedback gains are computed using the Riccati 
equation. The gain matrix that results is specific to the flight condition of the open loop model 
used to form the synthesis model. The next task is to define the feedback controller gains to 
meet the design requirements throughout the flight envelope. 

The longitudinal and lateral controller design examples found in Sections 7 and 8 respectively 
illustrate two different approaches to gain scheduling. The technique used for the longitudinal 
axis takes the Riccati solution for one middle-of-the-envelope flight condition and schedules it 
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based on flight condition parameters known to influence the airplane dynamics and control 
effectiveness (e.g., such as dynamic pressure and speed) in deriving a definition for gains that 
meets the requirements throughout the envelope. 

The gain scheduling technique used for the lateral controller example first solves the Riccati 

equation in determining the optimal gains at each flight condition. Whenever possible, the 

same set of criteria outputs and weighting matrices is used for all flight conditions. A check is 

made to make sure that the Riccati solution at each condition provides a closed-loop system 

that meets the design requirements. Gain schedules are then developed by plotting each gain ■ 

against a number of flight condition parameters by choosing the plot showing the greatest level 

of correlation and curve fitting the data. 

5.7 SENSOR SELECTION AND OUTPUT FILTERING 

The gain matrices produced in paragraph 5.6 include feedback gains from all of the synthesis 

model states to each of the control inputs. The next task is for sensor selection and output < 

filtering to provide the necessary feedback signals. State signals that are available directly 

from sensors are fed directly to the gain matrix. j 

Some states will be measurable over different frequency ranges using different sensors. For those 

states, complementary filters can be used to combine low frequency data from one source with ! 

high frequency data to another source. An example of this is altitude feedback where 

barometric altitude is used for low frequencies complemented by vertical acceleration data to i 

provide high frequency information. j 

Other states may not be directly measurable and require estimation. Although the design 
examples presented here do not include any such states, either full order Kalman estimation or 

i 

reduced order estimation could be used if the need arose. 

j 

5.8 CLOSED LOOP ANALYSIS 1 

j 

The chief objective of closed loop analysis is to verify that the design requirements have been 
met. In the event that the requirements are not met, closed-loop analysis should identify the 
problem areas to focus the design effort. Although different types of closed-loop analysis will 
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be used depending on the types of requirements in place, the analysis techniques described in 
paragraphs 5.8.1 - 5.8.6 will be sufficient to verify compliance with most design requirements. 

5.8.1 EIGENVALUES 


Computation of the closed-loop eigenvalues verifies compliance with mode frequency and 
damping ratio requirements. 

5.8.2 BROKEN LOOP FREQUENCY RESPONSES 

Gain and phase margins are determined via broken loop frequency analysis. The closed-loop 
system is broken, one loop at a time, at the input to each actuator and at each sensor output. 
Compliance with gain and phase margin requirements is verified by computing the margins for 
the resulting broken loop system. 

5.8.3 SINGULAR VALUE ANALYSIS 

In many cases the traditional robustness measures, gain and phase margins are not adequate 
because they are usually used on one loop at a time, thus overlooking the problem of 
simultaneous variations in different loops or, more generally, different system parameters. In 
their place, singular value analysis has been developed. Although singular value analysis is 
not used in the example designs presented in Sections 7 and 8, since it is not required by the 
design requirements (see Section 5), a brief description is included here for completeness. 

Singular value analysis is a multiple input/multiple output extension of classical single loop 
gain and phase margin analysis. Singular value analysis identifies the tendency of the system 
toward singularity (instability) for variations in more than one loop at a time rather than the 
one loop at a time structure for phase and gain analysis. An extension of singular value 
analysis, (i.e., structured singular value analysis) allows the designer to analyze the system for 
robustness to variations of specified structures. 

5.8.4 COVARIANCE RESPONSE 


Covariance response analysis predicts the RMS response of closed-loop system states and 
outputs to wind disturbance inputs. Wind turbulence is modelled as filtered white noise. The 
two most widely used turbulence models are Von Karman and Dryden filters. Each defines the 
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frequency content and magnitude of each component of six degree-of-freedom (three 
translational and three rotational velocities) air turbulence as a function of altitude, speed, 
and turbulence severity. 

After appending the appropriate wind model to the closed-loop airplane model, the 
covariance of the complete system to white noise is computed. The resulting response level data 
for airplane states, outputs, and control surface inputs are checked against the design 
requirements to verify that ride quality is within the required bounds and to control input 
activity levels. 

5.8.5 COMMAND FREQUENCY RESPONSE 

Closed loop frequency response analysis from pilot inputs to airplane states and control surfaces 
is used to check for compliance with performance requirements. Command response bandwidths 
must be at least as high as required without excessive control input activity (which might 
drive actuators into saturation) or high vehicle accelerations (exceeding structural or passenger 
comfort levels). 

5.8.6 TIME DOMAIN SIMULATION 

The final type of closed-loop analysis is time domain simulation. Often control system 
performance requirements are defined in terms of time histories. These requirements are 
verified via time domain simulation. In addition, all of the other analysis techniques are 
restricted to linear (or linearized) systems. Time domain analysis allows for inclusion of 
complex nonlinearities that cannot be properly treated using other analysis methods. 

! 

i 

I 

5.9 ITERATION TO SATISFY DESIGN REQUIREMENTS 

I 

An important feature of this, and any, design process is the iteration to satisfy all of the design 
requirements. The nature of control system design is making the trade-offs between control 
activity, performance level, and robustness to uncertainties and unmodelled dynamics. No | 

single pass design methodology would be able to sufficiently explore and optimize this trade- ! 
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5.10 SUMMARY OF DESIGN METHODOLOGY 


The strength of this controller design methodology is that it provides the power of LQG design 
in a way that is understandable from a classical frequency domain point of view. The intuition 
and experience of classical design is merged with the ability of LQG design to efficiently 
handle multiple input / multiple output design problems. Standard analysis tools have been 
employed to provide insight via open loop analysis. Historically proven frequency domain 
compensation techniques are used to form the criteria outputs completing the synthesis model. 
Once the synthesis model is complete, the Riccati equation is used to simultaneously solve for 
the complete set of feedback gains. After closing the loop, traditional analysis tools are used 
once again to verify compliance with design requirements. In addition, singular value and 
structured singular value analysis may be used to analyze multiple loop robustness. 
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6.0 DESIGN REQUIREMENTS 


The design requirements to be met by the controller consist of two sets: 1) the feedforward 
requirements based on ideal model specifications, passenger comfort and airplane dynamics 
constraints, and 2) the feedback requirements based on stability and performance robustness 
specification. These two distinct sets of requirements dictate the feedforward /feedback 
structure of the resulting controller. 

The design requirements for the longitudinal and lateral axes are presented separately in 
paragraphs 6.1 and 6.2, respectively. 

6.1 LONGITUDINAL AXIS DESIGN REQUIREMENTS 

A. Feedforward Loop 

1. Duplicate modes available in TECS 

a. speed hold (Mach, CAS) 

b. altitude hold 

c. glideslope capture and hold 

d. flare 

e. go-around 

f. flight path angle (FPA) 

2. Provide necessary command limiting to ensure 

a. passenger comfort 

b. proper energy distribution when throttles are at the limit 
B. Feedback Loop 

1. Robust Stability 

a . Gain margin of at least + 6dB, and phase margin of at least + 45 deg in all 
control and sensor loops. 

b. Sufficient rolloff in elevator loop at higher frequencies to avoid exciting 
unmodelled dynamics. 

c. Minimum phugoid damping ratio .55, minimum short period damping ratio .4. 
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2. Robust Performance 

a . Command loop bandwidth: 

• 06 - .12 rad/sec in cruise for both altitude and speed command loops 

• 1-1.2 rad/sec in landing for altitude command loop 

• .4 -.6 rad/sec in landing for speed command loop 

b. Control loop bandwidth: 

minimum required to achieve command loop bandwidth requirements. 

c. Gust response: 

same or better than TECS. 

It should be noted here, that in defining the requirements, TECS was used as a baseline, and the term 
"robust" should be interpreted to mean "throughout the flight envelope". 

6,2 LATERAL AXIS DESIGN REQUIREMENTS 

A. Performance Requirements 

1 . Zero Steady State Error 

a . Heading error when in heading mode. 

b. Localizer cross track error when following ILS localizer beam. 

c. Steady crosswind should not cause steady state errors in either mode. 

2. Zero Time Domain Overshoot 

a . Heading overshoot when in heading mode. 

b. Localizer beam overshoot when following ILS localizer beam 

(assuming capture is initialized with sufficient space to turn onto the beam without 
crossing the center line). 

c. Steady crosswind should not cause overshoot in either mode. 

3. Response Limitations 

a . Bank angle not to exceed 30 degrees during maneuvers. 

b. Lateral acceleration (nominally zero) not to exceed 0.05 g for maneuvers in still air. 

4. Bandwidth Requirements 

a. Maximize response bandwidths within limitations on overshoot, surface activity, 
eigenvalues, and disturbance responses. 
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B. Control Activity Requirements 


1 . Aileron Activity 

a . Position not to exceed +/- 15 deg 

b. Rate not to exceed +/- 30 deg/sec 

2. Rudder Activity 

a . Position not to exceed +/- 15 deg 

b. Rate not to exceed +/- 30 deg/sec 

C. Gust Response Requirements 

Gust reponse upper bounds given here are RMS values in response to a lateral Dryden gust with 
1 ft/sec RMS intensity. Units are feet and degrees where appropriate. 

1 . Heading and Sideslip Angular Responses 

a . Psi not more than 0.2 deg RMS 

b. Beta not more than 0.2 deg RMS 

2. Heading and Sideslip Angular Rate Responses j 

a . Psi not more than 0.2 deg/sec RMS 

b. Beta not more than 0.2 deg/sec RMS j 

I 

| 

3 . Aileron and Rudder Angular Position Reponses 

a . Aileron not more than 1 .0 deg RMS 

b. Rudder not more than 1 .0 deg RMS j 

4 . Aileron and Rudder Angular Rate Responses 

a . Aileron not more than 2.0 deg/sec RMS 

b. Rudder not more than 2.0 deg/sec RMS 

1 

i 

! 

5. Localizer Cross Track Response j 

a . Position not more than 5 ft 

b. Rate not more than 5 ft/sec j 

i 
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D. Robustness Requirements 


1 . Eigenvalues 

a . All eignvalues, 0.4 damping or better 

b. Dominant eignevalues, (lowest frequency), 0.6 damping or better 

2. Stability Margins 

a . Aileron input: simultaneous +/- 6dB and +/- 45 deg 

b. Rudder input: simultaneous +/- 6dB and +/- 45 deg 

c. All sensor inputs: simultaneous +/- 6dB and +/- 45 deg. 


7.0 LONGITUDINAL AXIS DESIGN 


The longitudinal controller structure reflects the design requirements discussed in paragraph 6.1 and 
consists of the feedforward and feedback loops. 

The feedforward controller receives pilot commands through the Mode Control Panel (MCP) and 
generates altitude, altitude rate, speed and speed rate commands, which the feedback regulator 
loop must track. The decision to have both position and rate commands was made in order to follow 
TECS structure as close as possible, and also to be compatible with 4D mode (time and space) if one is 
added in the future. 

The feedback regulator receives both feedforward commands and airplane sensor signals from which 
it generates elevator and throttle commands. Figure 5 shows the general structure of the control 
system as well as all the major interfaces. The pilot, via MCP, selects the mode he wants the 
airplane to fly by entering the desired altitude and speed targets. At any given time, only two 
modes can be engaged: one speed mode and one altitude mode (i.e., each mode generates either speed 
or altitude commands). Hence, the feedback regulator has a common interface with all feedforward 
modes. 


7.1 FEEDFORWARD MODES 

Each feedforward mode shown in Figure 5 will be described in detail in the following paragraphs. 
7.1.1 ALTITUDE HOLD MODE 


The altitude hold mode is a direct engage mode, (i.e., when pilot engages this mode any other mode 
that controls altitude is automatically disengaged) consisting of a second order filter which limits 
the first and second derivatives of altitude command (i.e., altitude rate and acceleration). The 
bandwidth of the filter is selected to assure proper altitude command response. The filter receives 
an altitude target from MCP and generates altitude and altitude rate commands for the feedback 
regulator as shown in Figure 6. The transfer functions from altitude target to filter outputs are: 


F(Hc,HmCP) = 


sW n 2 

s 2 + 2^W n s + W n 2 


and 


( 1 ) 
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MCP LOGIC FLAGS 


- 

ALT 

HOLD 



- 

GS 

CAP & 
HOLD 




FPA 



- 

GA 



- 

FLARE 



- 

SPD ■ 
HOLD 



ELEVATOR COMMAND 


THROTTLE COMMAND 


AIRPLANE SENSORS 


Figure 5. Controller Structure 



HUM = ± (IHDOTI+IO) W n =.15 

H LIM = O.OSg = 1.6 fps* ; = 1 

[> IF (.NOT. ALTH) 1 1 = H LIM TECS 

[5> IF (.NOT. ALTH) 12 = H 


Figure 6. Altitude Hold Mode 
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( 2 ) 


H c _ W n 2 
HmCP s 2 + 2£W n s + W n 2 

where 

^ = 1 


W n = .15 


The altitude rate limit (Hum) taken from TECS is rather loose. Hum = ± ( I Hum I + 10) , hence 
the H signal will rarely be on the limit. (HfJM ) is much stricter: Hum = ± *05 g and is introduced 
to ensure passenger comfort. The value of this limit was selected so as to never command a total of 
more than .lg during simultaneous altitude and speed command changes, since speed command second 
derivative is also limited to .05g. 

When the airplane is not in the altitude hold mode, the H integrator (II) should be initialized to 
current altitude rate, and the H integrator (12) to current altitude. 

It should be noted that the natural frequency of the filter was selected to be slightly higher than 
the required altitude command bandwidth of .06 to .12 rad/sec (W n = .15). W n was adjusted to ensure 
an overshoot free linear response of the total control system to pilot commands. Since in nonlinear 
situations the filter will be on the H limit most of the time, the value of W n is not of great 
importance. 

7.1.2 GLIDESLOPE CAPTURE AND HOLD 

The glideslope mode shown in Figure 7a is a direct arm/automatic engage mode (i.e., when the pilot 
receives a valid glideslope signal he arms the glideslope mode). Glideslope capture is 
automatically engaged and disengages any previous mode when a certain criteria is satisfied. 

To understand how the glideslope mode works, one must have a clear picture of glideslope geometry 
(i.e., the position of the airplane with respect to glideslope beam and the variable definitions 
associated with it). The glideslope geometry is shown in Figure 7b. 

The valid glideslope signal received by the onboard receiver is the glideslope error (GSE) or the 
angular deviation of the airplane’s flight path from glideslope beam centerline. Based on 
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DIRECT ARM / AUTOMATIC ENGAGE MODE 



GLIDESLOPE CAPTURE CRITERIA: 

IF ((HE/15 + (H c - H)) HE/15 .LT. 0) FAGLD = 1 


GLIDESLOPE HOLD CRITERIA: 

IF ( IHEI .LT. 10 .AND. | H c - H I .LT. 3) FAGLDH = 1 


Figure 7a. Glideslope Capture and Hold Mode 



Figure 7b. Glideslope Geometry 
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glideslope geometry, the GSE angle is converted into altitude error (H e ) and the valid altitude 
command is then computed: 

Hc = HRAD + H e (3) 

The H command is computed to avoid switching inside the feedback loop, (since GSE is an error 
signal) and maintaining the uniformity of the feedback/ feedforward interface. Its computation is 
based on the assumption that the glideslope angle is known and usually equal to 3 degrees. 

Computation of H e and H are presented in Appendix A. 

The glideslope capture criterion is based on the feedback regulator structure, namely, the altitude 
integrator computation, shown in Figure 7c. 



Figure 7c. Altitude Integrator Command Computation 


The airplane is on the overshoot-free glideslope capture envelope when H e = 0 , or 

C C 


.1 H e + 1.5 H e = 0 

or (4) 

l/15H e + He = 0 

When (4) is zero, the expression H e /15 • (H e /15 + H e ) changes sign. This criterion is used to set 
glideslope engage flag (FAGLD). 
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The feedback regulator has two separate sets of gains: one for cruise and one for landing. The switch 
between the two occurs when the airplane is on the glideslope beam. The successful capture of the 
beam is determined by the following criterion: I H e 1 <10 ft and IH e I <3ft/sec. When this 
criterion is satisfied, the glideslope hold flag (FAGLDH) is set, and the landing gains replace the 
cruise ones. 

7.1.3 FLARE MODE 

The flare mode is an automatic engage mode and is engaged at 45 feet above the runway. The mode- 
generates an inertial altitude path to land the airplane 1200 feet from the flare engagement point as 
shown in Figure 8. 

An extensive study of autoland flare control has been done by A. Lambregts et al for the NASA TSRV 
airplane [7]. The study concluded that inertial path tracking is a better way to accomplish low 
touchdown dispersion. This recommendation has been followed in this design. 

An inertial path must satisfy four boundary constraints: at flare engagement, the altitude command 
must equal 45 feet and the altitude rate command must equal the current descent rate of the airplane, 
and at touchdown, the altitude command must be zero and altitude rate command must equal -2.5 fps. 
The algebraic function selected to satisfy these constraints was a hyperbola, as depicted in Figure 8. 
The computation of the polynomial coefficients A, B, C and D can be found in Appendix A. 

7.1.4 GO-AROUND MODE 

The go-around mode is a direct engage mode, used in case of an aborted landing. The airplane is 
commanded to climb out at 10 deg FPA. As shown in Figure 9, FPA command (y c ) is converted into 
H c . as follows: 

He = Yc • Vi/57.3 (5) 

Hc = jH c (6) 

The integrator in (6) should be initialized to current altitude when not in go-around mode. 
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H c (x) = Ax 3 + Bx 2 + Cx + D 

V V V 

H c (x) = 3A-5J-!- x 2 + 2By-i- x + C-^ 


TD 


TD 


TD 


BOUNDARY CONDITIONS: 


H c (o) = 45 ft. 


H (°) = H of airplane 


H c( x td) = 0 


»c( x td) = -2- 5 


Figure 8. Flare Mode 



[i> IF ( .NOT. FAGA ) H = H 


Figure 9. Go - Around Mode 


1200 ft 


c 
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7.1.5 FLIGHT PATH ANGLE MODE 


The Flight Path Angle (FPA) mode is a direct engage mode. The pilot enters the desired FPA he 
wants the airplane to follow through the MCP. Commanded FPA is then converted to He and He, as 
described in the previous sections, and then processed through a second order filter as shown in 
Figure 10. This filter is similar to the one used for the altitude mode, except He (computed from y c ) 
comes through a feedforward path to improve filter command tracking. The FPA mode was designed 
to duplicate TECS capability, but its intended function can be successfully accomplished by using 
altitude hold mode. 

7.1.6 SPEED HOLD MODE 

The Speed Hold mode is a direct engage mode. The Pilot enters the desired speed target, MACH or 
CAS, into the MCP. The command is then converted into TAS command (see Appendix A) as shown 
in Figure 11. The TAS speed target is then passed through a second order filter to generate speed and 
acceleration commands, as required by feedback regulator interface. The filter dynamics are 
identical to the altitude hold filter dynamics. Thus, when in the linear region, both filters provide 
coordinated commands to feedback regulator. 

The speed command filter has V and V limiters (Figure 11). TheV limiter is computed to provide 
vertical acceleration limiting of .05 g ^ yj^S j to ensure passenger comfort. TheV limiter is used 

to limit forward acceleration to .lg (3.2 fps) when the throttles are in the linear region. When the 
throttle is on either limit, the feedback regulator is in the speed-on-elevator (SOE) configuration, 
and V limit is computed to ensure proper energy distribution between potential and kinetic energies 
of the airplane. When throttles are at the forward limit, the upper value of Vjjm is computed to 
ensure that a minimum rate of climb (10 fps) is maintained. When throttle is at the aft limit, the 
lower value of Vlim * s computed to allow the airplane to level off, but not climb in descent. The 
formulas for both limiters are explained in Appendix A. 

When not in the speed hold mode, both integrators (II and 12) should be initialized to current true 
acceleration and VTAS, respectively. 
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H LIM a ± ( I H | +10) 

HUM a O.OSg s 1.6 Ips 2 
[> IF (.NOT. FAFPA) I 1 = H 

[£> IF (.NOT. FAFPA) I 2 a H 


W n a .15 

C ■ 1 


I> 


IF (.NOT. FAFPA) I 3 = H 




I 


Figure 10. Flight Path Angle Mode j 
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Figure 11. Speed Hold (Mach or CAS) Mode ' 
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7.2 FEEDBACK REGULATOR DESIGN 


Reference [3] describes in detail the methodology used to design the feedback regulator. Here, this 
methodology is applied to design a controller to satisfy the feedback requirements. 

It is clear that there are two distinct sets of requirements for the cruise and landing portions of flight 
envelope (i.e., a more sluggish airplane is desired in cruise, whereas in landing a tight tracking of 
glideslope beam and flare path is required). Hence, the decision was made to design two sets of 
feedback gains: one for cruise and one for landing. The structure consists of an integral regulator plus 
a complimentary filter for true airspeed (VTAS) and acceleration. The altitude (H) and altitude 
rate (H) signals are already synthesized by the onboard IRU. 

The integral regulator portion of the design will be presented first, assuming V and V are available, 
followed by a description of the complementary filter. 

The standard transport airplane longitudinal dynamics model consists of four states: u, a, Q, 0 and 

two control effectors, 8 e and 5 th* To complete the model two more states were added: H and EPR 

(an output of first order engine model as shown in Figure 12). Since a is a very noisy signal, it was 

VTAS • 

replaced by H state, using a similarity transformation: a = 0 - ^ ^ H . All flight conditions 
used for linear controller design and analysis are listed in Appendix B. 

The following paragraphs describe the design of the cruise and landing regulators. 

7.2.1 CRUISE REGULATOR DESIGN 

For cruise design, a heavy weight cruise condition (#159) was selected. It is listed in Appendix A. 
The design then proceeded in the following steps: 

1 . Open loop analysis 

2. Output criteria creation 

3. Diagonal weightings, Q and R, selection 

4. Closed loop analysis 

Each step, or the combination of several of them, were iterated many times before a satisfactory 
solution was achieved. Results of each step are now briefly described. 


41 


7.2.1.1 Open Loop Analysis 


The results of open loop analysis can be found in Appendix A. There were no surprises. Both the 
elevator and throttles have sufficient gains at d.c. to provide steady state control of either altitude 
or speed, as is clear from SISO frequency responses. The elevator is the best effector to control both 
the phugoid and short period modes. The throttles best control the energy mode. There is right half 
plane transmission zero at 7 rad /sec, which is outside of the frequency range of interest, and hence is 
of no concern. There are no observability problems, since all the states are available. 

7.2.1.2 Output Criteria Creation 

The complete synthesis model (airplane plus output criteria) is shown in Figure 12. Two criteria 
outputs were created, H_CRIT and V_CRIT (see Figure 12) because there are two control effectors 
(elevator and throttles). Furthermore, since the synthesis model is square, the transmission zeros of 
the total system are the ones of the airplane plus the ones created by criteria outputs. 

H_CRIT output is a frequency weighted combination of altitude (H), altitude rate (H) and vertical 
acceleration (H) and has the following expression: 

Kt 

H_CRIT = y [ .1 (-H c + H) - He + H ] + Kpi (-H c + H) - Kp2 H + H (7) 

Since 


H = sH 
H = s 2 H 


then 

H CR1T = ■ K1S2 1 KlS * 1Kl He , * K1S2 ± K2S + - K 3 h 
s s 


where 

Ki =Kpi = 3.2 
K 2 =Ki- K P2 = 5 

K p2 = 5 
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Figure 12. Cruise Controller Synthesis Model 










K 3 = .1Ki = 1 


Therefore, H_CRIT adds three transmission zeros to the synthesis model, which are the roots of the 
polynomial s^ + 3.2s^ +5s + 1 = (s + .232) (s^ + 2.97s + 4.31). 


The real transmission zero will attract the altitude mode, and the complex pair of transmission 
zeros will attract the short period of the aircraft. 


The V_CRIT output is a combination of airspeed (V) and acceleration (V). The expression for 
V.CRIT is: 

Ki . . . 

V_CRIT = -j [(.1(V C - V) + V c - V) ] - KpV + V c - V or 


s^ + KjS + .lKi s^ + (Kp + .1 Kf)s + .lKi 

s c s 


(9) 


where 

V =sV 
K p = .3375 
Ki = .625 


V_CRIT adds a pair of complex zeros to the synthesis model to attract the phugoid mode of the 
airplane. The zeros are the roots of the polynomial s^ + .4s + .0625. All the zeros of the synthesis 
model are listed in Table 5. Once the criteria variables have been selected, diagonal elements of the 
Q & R matrices in the cost function J must be determined, based on desired crossovers in broken loop 
responses, [3]. J assumes the following form: 
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CRIT V_CRIT 8 P 5 


TH, 


QO' 

OR 


fH_ CRIT' 
V CRIT 
8e c 
8. 


TH 


c J 


( 10 ) 


The values of Q & R matrices and resulting feedback regulator solution are shown in Table 6. 


7.2.13 Gain Scheduling 

Once the nominal cruise regulator had been designed, it was tested for all cruise and glideslope 
capture conditions (#1 - 79, 113-160). The results showed that gain-scheduling is necessary to 
achieve uniform response throughout the flight envelope. 


The following schedules were used: 


8 


• nominal 8 e 

0 


250- W 


TH c 110,000 KTHR 


nominal 8 


TH„ 


KH = 1.5 • nominal KH 


where 


Q - dynamic pressure 
W - airplane’s weight 


KTHR - thrust to throttle handle gain ratio 


KH - altitude rate feedback gain on the H input to the controller 
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Table 5a. Transmission Zeroes of Cruise Controller Synthesis Model 


REAL 

IMAG 

DAMPING 

FREQ 

6.964 

0 

1 

6.964 

-.2319 

0 

1 

2.319 

- .2 

- .15 

.8 

.25 

- .2 

.15 

.8 

.25 

- 1.484 

-1.452 

.7147 

2.076 

- 1.484 

1.452 

.7147 

2.076 

- 6.888 

0 

1 

6.888 


Table 5b. Cost Function 


J = 


oo 

f TH CRITV CRIT 5„ 8 TU 1 [ Q Ol 
0 L c TH c J L O R J 


H_ CRIT 
V_ CRIT 

8 -C 

S TH. 


dt 


Table 6a. Weighting Matrices for Cruise Controller 

H_CR1T V_CR1T 


q = r 20 ° i 

L o I.6E-4J 

^TH 

LO 0.6j 


Table 6b. Feedback Gain Matrix for Cruise Controller 


HDOT 


THETA 


H 


EPR 


r»..i 


.1787 

6.01 E-02 

.6223 

.5769 

1.53E-03 

-.4178 

.2871 

1.9E-2 


1 

8 th 

c 


-.7494 

-.5444 

.1574 

-2.865 

-.1143 

-6.067 

2.981 

-9.2E-2 


1 

i 
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7.2.1.4 Speed on Elevator (SOE) Configuration 


An important issue concerning the feedback regulator design is the throttle limiting cases. When the 
throttles are either at the forward or aft limit, the elevator will control airspeed, and, as described 
earlier, speed command processor will properly distribute aircraft’s total energy. Therefore, when 
the throttles reach the limit, the altitude integrator and altitude gains in the elevator command 
computation are set to zero (Table 6). This design results in poor phugoid damping for the new closed 
loop system. To improve the damping, both speed integrator gains are increased 2 times, and the 
total elevator command gain by 2.5 times. 

Figure 13 shows the general structure of the cruise controller. It is similar to TECS structure [4], 
except the cross-couplings are at the outputs of the integrators, rather than at the inputs as in TECS. 

7.2.2 LANDING REGULATOR DESIGN 

The landing regulator was designed to satisfy glideslope hold and flare requirements and hence, it 
has much faster command responses accompanied by much higher feedback gains than does the 
cruise controller. 

The nominal design of the landing controller was done for flight condition #102. This is a landing 
condition with flaps set at 30 degrees and landing gear down. 

The open loop analysis (Appendix A) shows a right half plane transmission zero at 3 rad/ sec, inside 
the frequency bandwidth trying to be controlled. This zero is one of the pair of so-called "percussion 
zeros**,- usually nea rly symme tr ic about imaginary axi s . -The non -minimum ph a se ze rods due to the 
fact that airplane's center of percussion is forward of the c.g. If the N z and H sensors are moved 
forward of the center of percussion, then the zeros will become minimum phase. This is done 
mathematically using the formulae: 

N z — N ZC g + L/G Q 

H = Hcg + L • Q (ID 

H = Hcg + L • 0 

where L is the distance forward of the center of gravity. 
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Figure 13. General Structure of Cruise Feedback Controller 
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Figure 14. Landing Controller - Synthesis Model 



















Figure 14 contains the synthesis model for the landing controller with redefined H, H, and H as in 
(11) ( i.e., the H and N z sensors have been moved 6 feet forward of the airplane’s c.g.). 

The transmission zeros created by criteria outputs, the Q and R weighting matrices and the regulator 
feedback solution for landing controller can be found in Table 7. The mathematically redefined 
percussion zeros of the airplane are now at 13.09 rad/sec and -14.61 rad/sec, outside of the control 
bandwidth. 


7.2.2.1 Gain Scheduling 


Once the nominal design for the landing controller had been completed, it was tested on the 
remainder of landing conditions (# 1-112). The results were found to be inadequate; therefore, the 
following gain scheduling scheme was developed. For each landing condition, a regulator solution 
was obtained for the cost function used to design the nominal controller. Thus, a total of 16 feedback 
gains had been derived for each of the 112 conditions. Each of the gains was then plotted versus a 
flight condition parameter (e.g., such as H, a, Q, KTHR, etc.). 


Figures 15 and 16 present an example of such plots. In Figure 15, Kj-j (altitude to throttle command 
gain) is plotted versus Q (dynamic pressure). For each flap setting there is a clear hyperbolic 
dependence of Kh on Q. Figure 16 shows a plot of K\j (airspeed to elevator command gain) versus Q. 
A hyperbolic relationship between K\j and Q is obvious and is independent of flap setting. The 
expression of the form: 


K = 


A + BX 
1+CX 


( 12 ) 


was used to find a curve fit for each gain. In (12), A, B, and C are either constants or functions of flaps 
and X is an independent variable (e.g., Q, a, H, etc.). Appendix A contains complete information for 
each gain. 

7.2.2.2 Airspeed Complementary Filter 

The airspeed filter (Figure 14) uses airdata and inertial signals to generate accurate airspeed and 
acceleration signals over a large frequency range. As shown in Figure 14, the filter's time constant x 
is scheduled as a function of altitude since tight tracking of airspeed is not required at higher 
altitudes. 
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Table 7a. Transmission Zeroes of Landing Controller Synthesis Model 



REAL 

IMAG 

DAMPING 

FREQ 

13.09 

0.0000 

-1.000 

13.09 

-0.8900 

0.9818 

0.6716 

1.325 

-0.8900 

-0.9818 

0.6716 

1.325 

-1.200 

0.9000 

0.8000 

1.500 

-1.200 

-0.9000 

0.8000 

1.500 

-1.777 

0.0000 

1.000 

1.777 

-14.61 

0.0000 

1.000 

14.61 


Table 7b. Weighting Matrices for Landing Controller 


H_CRIT V_CRIT 

r 1.5 o i 


R = 


5 e 5 TH 

10 0 
0 10 


Ci 


Table 7c. Feedback Gain Matrix for Landing Controller 




V 

HDOT 

Q 

THETA 

H 

EPR 

1 

V 

>h 

rv. 1 


0.6035 

2.733 

3.491 

7.425 

2.173 

8.624 

0.04 

1.2 

c 


-1.9 

-0.3028 

-0.2504 

-0.5707 

-0.4037 

-17.85 

0.3876 

-0.1263 
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THRUST GAINS SCHEDULE 




Figure 16. Landing Controller - Elevator Gains Schedule 







73 DIGITAL IMPLEMENTATION 


The feedback controller (with the regulator solution and complementary filter solutions combined) 
has 3 states (complementary filter state, altitude and speed integrator states), 7 sensor, 4 command 
inputs and 2 outputs (elevator and throttle commands) and the following general form: 


x = Ax + By 
u = Cx + Dy 


(13) 


where 


A = 


-J- 0 0 


(4- 


1J 0 0 

0 0 0 


(13a) 


0000 ^ 10 0 0 0 


B ^0 0 0 0 -^- 1. 1 0 1 


.1.5 0 0 0.1 0 0 -.1 0 - 1.5 


(13b) 


nc„ k. 


c = 


K, 


8 fi V 5 

e 8 . 


K, 


K, 


TH 


H. 


K, 


H. 


TH 


THJ 


(13c) 


D = 


pK h 
H * e 



EPR S 

5 e 

k h 

H «e 

K v 

5 e 

0 

0 

0 

0 

0 

K *s 
L TH 

X. 

k \ h 

k epr 8 

TH 

KH 8 

TH 

K v 

6 th 

0 

0 

0 

0 

0 


(13d) 
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u = (H, Q, 0, EPR, H, VTAS, Vi, V c , Ho V c< He ) T 


y=(8e t -5 THc )T 


where the K’s in the C and D matrices are distinct for cruise and landing solutions. 
The feedback controller was discretized using the following formula: 


x — p AAT v _l 
X k+1 e X k + 


I 

V 0 


! A<7 da 


Bu. = A , x. 4- B , u, 
k d k d k 


yk+l=Cxk + Duk (14) 

where 

AT = sampling period 

xk = discretized controller state vector 

Since x is a function of altitude in 13a and 13b, closed form expressions had to be obtained for Ad and Bd in 
(14). The results are presented in Appendix A. 


The discretized feedback controller was implemented digitally using delta-coordinates concept (Figure 18). 
The 's’ and T/s’ meaning differentiation and integration, respectively, are an abuse of notation in this case. 
Rather than computing actual derivatives of the command and sensor signals, the difference between present 
and previous values of each signal is found. At the output, the integrator is initialized to current elevator 
and throttle positions, and the controller outputs are then added to the previous values of the integrator. 



COMMAND 

COMMAND 


Figure 18. Digital Implementation - Delta Coordinates 

The delta concept offers significant advantages over the full signal implementation: 

a . During controller switching (e.g., between cruise and landing controllers) all states of the new 
controller should be initialized to 0. 
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b. Full scale implementation may lead to incorrectly computed gains, which are scheduled as a 
function of quickly changing flight parameters (e.g., a). This problem is avoided with delta 
concept [8]. 

c. The issues that arise during throttle limiting can be easily dealt with as described in the 
following. 

| 

I 

Figure 19 shows in greater detail the delta implementation of the feedback controller. When ( 

FAGLDH flag is set, the switch between cruise and landing controllers takes place. At this point, 
all landing controller states are initialized to 0. Both controllers generate an incremental throttle 
command, § . When the throttle lever reaches the forward limit (FWDLIM) and 8 jjq is 

positive, it should also be set 0. If the throttle lever is at the aft limit ( AFTLIM) and 8 is 
negative, it should be set 0. By setting 8 to 0, throttle integrator saturation is avoided. Also, 
the sign of 8 determines when to take the throttle off the limit. The advantages of this 

c j 

strategy are described in paragraph 7.4.2. I 

| 

Some of the feedforward loop modes have second order filters, which were also discretized using j 

expression 14. It should be noted that the feedforward loop modes were not implemented using a 

delta coordinates concept. The Aj, Btf matrices for the discretized feedforward loop are given in 

Appendix 7. ! 

i 

7.4 RESULTS j 

i 

i 

< 

The complete analysis and testing of the total design consisted of two parts: 1) linear analysis of the ] 

feedback controller, and 2) nonlinear testing of several feedforward modes, (i.e., feedforward and 

feedback loops combined). \ 

7.4.1 LINEAR ANALYSIS 

I 

Linear analysis of the feedback regulator included the following items: ^ 

a . Open loop eigenvalues j 

i 

b. Closed loop eigenvalues j 

l 

i 

i 
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Figure 19. Throttle Limiting - Controller Switching 



Appendix A contains the results of the linear analysis of the cruise controller over a substantial 
portion of the flight envelope. 

Plots 1-4 show open and dosed loop eigenvalues for 48 (#113 - 160) flight conditions. It is clear that 
significant improvement in eigenvalue damping has been achieved with feedback. Both phugoid 
and short period modes damping ratios satisfy the closed loop damping requirements. 

Plot 5 shows gain margins plotted versus flight condition for conditions # 113-160 for all control and 
sensor loops. Plot 6 shows phase margins for the same conditions. It is clear from both plots that 
gain and phase margin requirements have been satisfied. 

Plots 7 and 8 show crossover frequencies for elevator and throttle control loops. The maximum 
elevator loop crossover frequency is at 4.1 rad/sec, which is well below 15 rad/sec (the bandwidth of 
the elevator actuator). The throttle loop crossover frequency remains constant throughout the flight 
envelope at around .35 rad /sec, which is, again, well below 1 rad/ sec (commonly accepted to be the 
engine bandwidth). 

Plots 9-12 show closed loop frequency responses of the airplane’s altitude and speed to altitude and 
speed commands for conditions #113 - 128. There is little variation in responses as flight conditions 
change. In fact, this sample of 16 conditions is representative of the rest of the flight envelope. The 
H/H c (Plot 9) response has a bandwidth of .08 rad/sec, V/V c (Plot 12) has a bandwidth of .1 
rad/ sec. Both are within the command response bandwidth requirements for the cruise controller. 
Plots 10 and 11 show the crosscoupling effects of altitude and speed commands. H c has little effect 
on speed response (Plot 10). V c effects the airplane's altitude significantly in the frequency range 
between .01 and .5 rad/sec. This result is to be expected because of the low bandwidth feedback 
controller in cruise. This kind of altitude response is also wanted by the pilots, who would rather 
have the airplane drop 100 feet, than see throttles move 1 degree. 

Plots 13 - 19 show closed loop aircraft covariance responses to 1 fps longitudinal and vertical Dryden 
turbulence. Plots 13, 14 and 15 show airspeed, altitude and N z responses plotted versus flight 
condition. All show very small variations around the mean. Plots 16 and 17 show elevator and 
throttle covariance responses and plots 18 and 19 show elevator and throttle rate responses to 
Dryden turbulence plotted versus flight condition. It is clear from the plots, that there is little 
control effector activity in turbulence during cruise. 
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7.4.2 NONLINEAR TESTING 


The term 'nonlinear' is an exaggeration, since only the nonlinear controller was implemented, the 
airplane model remained linear. The nonlinear controller implementation included discretized 
linear controller with all the limiting and switching which takes place in the feedforward loop. 
This implementation was done to illustrate the nonlinear capabilities of the controller. The results 
of the 'nonlinear' simulation are presented in Appendix A. 

Plots 20 - 24 present the results of a simulated glideslope capture and hold. Plot 20 shows the 
altitude and glideslope profiles. The discontinuity in glideslope profile indicates the time when 
glideslope capture mode is engaged. Proper switching time results in overshoot free capture of 
glideslope beam. Once the beam is captured, a switch occurs from cruise to landing controllers. 

Plots 21 - 24 show the rest of the airplane's variables to better demonstrate the airplane's 
performance during capture. 

Plots 25 - 29 demonstrate the airplane's performance in speed hold and altitude hold modes. The 
airplane is trimmed at 35,000 feet when the pilot dials in a simultaneous altitude change command 
of 3,000 feet and speed increase command of 30 fps, conventionally known as flight level change 
maneuver. 

Plot 25 shows the altitude profile and points out throttle limiting times. The simplicity of the logic 
associated with throttle limiting is obvious. 

Plot 26 presents the same profile, except the logic events are shown that would normally take place 
iruthe conventional control-system , consis t ing of sepa r a te autopilot and autothrottle boxes. Once a 
sufficiently large altitude change is requested, the autothrottle switches into MAX EPR mode to 
control the engines to maximum EPR for given flight condition. The maximum EPR command is 
retrieved from the data base. The autopilot switches to speed on elevator (SOE) mode. While the 
aircraft is climbing, the autopilot computes the top of climb (T/C) point where it will switch to a 
path capture mode, which uses a nonlinear capture controller. At this time, the autothrottle will 
switch back to speed on throttle mode. Once the desired altitude target is captured, the autopilot 
will switch to a path hold mode, a linear altitude hold control law. The operation of the 
conventional control system is presented here to demonstrate the advantages of the integrated 
design. 
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Plot 27 shows airspeed, altitude rate and throttle position profiles for the flight level change 
maneuver. These signals are plotted together to better illustrate the limiting that takes place in 
the speed command processor when throttles reach the forward limit. Throttle limiting occurs after 
10 sec into the maneuver. As shown, airspeed response levels off at 25 fps above trim value and 
altitude rate levels off at 10 fps, while throttles are at the limit. This behavior is a result of V 
limiter in the speed command processor (i.e., when throttles are at the limit the speed command 
processor limits the speed command to maintain a minimum rate of climb of 10 fps). This is the 
reason why airspeed levels off at 25 fps, rather than reaching the target of 30 fps right away. 

When the throttles go off the limit, the airplane starts to pitch over and the remainder of the speed 
target and altitude target are captured. In addition, the ripple created when the throttles initially 
go off the limit in throttles position profile is due to the fact that the excess energy released by the 
airplane pitching over is used to capture the remainder of the airspeed target. 

Plot 28 shows the change in N z level of the aircraft during the flight level change maneuver. It 
doesn't exceed .05 g, which was achieved by the H limiter in the altitude command processor. 

Plot 29 shows the elevator position profile for the same maneuver. 
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8.0 LATERAL AXIS DESIGN 


This section presents an airplane autopilot design process, illustrated by an example design for a 
transport airplane, using an integral LQG design technique. Fourteen linear models trimmed at 
different flight conditions were used during design and analysis phases (see Section 4 for open loop 
airplane description). 

The objective was to design a single autopilot control system to provide both lateral axis stability 
augmentation and aircraft directional control across the flight envelope. The autopilot provides 
heading and ground track heading hold for cruise, and localizer beam capture and hold for landing 
approach. 

The approach was to use ailerons and rudder to independently control heading and sideslip. The 
feedback controller was designed to provide stability augmentation and sufficient command response 
bandwidth to meet the performance requirements. A separate feedforward controller was designed 
to filter pilot inputs to achieve desired transient responses. An additional outer loop controller 
generates heading commands for localizer capture and hold. 

8.1 HEADING CONTROLLER DESIGN 

The heading controller was designed using an integral LQG model following technique. (Figure 20 is 
a block diagram of the controller structure.) The feedforward and feedback controllers were 
designed separately. Feedforward design was driven by performance requirements while feedback 
design was dominated by stability and robustness issues. The following paragraphs outline the 
design q^rocess referring to the specific example for illustration. 

8.1.1 REGULATED VARIABLES 

For each independent control input a regulated output is chosen. For this example the independent 
inputs were aileron and rudder while the regulated variables were heading and sideslip. Likewise, 
for each regulated output a pilot input is designated and an ideal model defining the desired 
response of the regulated variable to pilot inputs is selected. 

For this design the input for heading command was a compass direction entered via the autopilot 
control panel. Sideslip command was given by the rudder pedals (even though it is rare that a 
sideslip command would be input during autopilot flight). The ideal model for heading response 


61 



Figure 20. Closed Loop Block Diagram of Integral Model Following Structure 
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8.13 GAINS AND GAIN SCHEDULES 


The steps given in paragraph 8.1.2 are part of an iterative process that involves solving the Riccati 
equation for the feedback gain matrix followed by closed loop analysis to check for compliance with 
the design requirements. The objective is to find a single set of criterion outputs and weighting 
matrix values that yield adequate closed loop characteristics throughout the flight envelope. 

Failing this, the flight envelope may have to be partitioned into sub-regions and a separate design 
completed for each. The example presented here yielded a single set that proved satisfactory 
throughout the envelope. 

Having defined the criterion outputs and weighting matrix values, the next step was computation of 
the feedback gains at each of the fourteen flight conditions presented with this problem. The results 
(see Figure 23) show that some of the gains remained relatively constant for the various flight 
conditions while others exhibited large variations. Since the airplane flies through a continuum of 
conditions, a method of defining the gains to be used at any flight condition is necessary. 

The task of generating continuous functions defining the gains is called gain scheduling. Gain schedules 
define the gain values as a function of measurable flight condition parameters (e.g., dynamic pressure, 
speed, altitude, flap angle setting, etc.). Each gain is scheduled separately by plotting the design 
values against the various flight condition parameters. The flight condition parameter exhibiting the 
greatest correlation is used for scheduling a particular gain. The schedule function is derived by curve 
fitting the plot of design values versus flight condition parameter. For the heading controller design 
example, some of the gains were set to constant values while others were scheduled against calibrated 
airspeed, mach number, or flap setting (see Figure 24). Note that Figures 23 and 24 the integral Beta to 
aileron gain Itas been^eroed^This was donejo avoid instability in the event of a loss of rudder. 
Zeroing this gain had little effect on closed loop performance and robustness. 

Results are shown here for fourteen flight conditions. Plots 30A and 30B found in Appendix B along 
with the accompanying gain schedules and controller block diagram define the controller for the 
complete design suite of 160 flight conditions. 

8.1.4 BETA COMPLEMENTARY FILTER 

The final step in the design process is estimation of any feedback states that are not directly 
measurable with the required fidelity. If needed, a full state Kalman estimator can be designed. 
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was a 0.4 rad /sec lag including a rate limiter to limit bank angle to 30 degrees. The sideslip ideal 
model was a 0.5 rad/sec lag. Error signals used for feedback control are defined by taking the 
difference between the output of each ideal model and its corresponding regulated variable. 

(Figure 21 is a block diagram of the feedforward controller.) 

8.1.2 CRITERION OUTPUTS AND WEIGHTING MATRICES 

The next step is to form criterion outputs for use during LQG synthesis. One criterion output is defined 
for each regulated variable. 

The first step in forming the criterion outputs is to provide high gain at low frequency between the 
control inputs and the criterion outputs. The result is large penalties on steady state errors leading to 
good steady state tracking. Often, as was the case between aileron and heading, the plant itself 
provides sufficient low frequency gain. In other instances, as with sideslip, an integrator is 
introduced to provide infinite gain at steady state. 

The second step in forming the criterion outputs is to examine the transmission zeros of the synthesis 
model. Because the synthesis model is square (same number of control inputs as criterion outputs) its 
transmission zeros can be computed. An important feature of this design technique is the asymptotic 
tendency of the closed loop eigenvalues to migrate toward the transmission zeros. The plant itself 
will have transmission zeros over which the designer has no control (other than chosing different 
inputs and outputs). In addition to the natural zeros, the designer can modify the synthesis model by 
adding other output signals to the criteria outputs to create additional zeros. In the example 
presented here, a real zero was added to the heading criteria (s = -3) to attract the heading state. A 
complex pair of zeros (Wn = 2, Zeta = .8) was added to the sideslip output to attract the dutch roll 
mode. 

Once the criterion outputs are defined, the final step is selection of the diagonal Q and R weighting 
matrices for LQG synthesis of the feedback gain matrix. Q and R have the same dimensions since 
there are the same number of control inputs and criterion outputs. In this example both were [2x2], 

As with any LQG technique, the weightings are chosen to make the trade off between control and 
output activity. To aid in choosing the weightings the frequency responses from control inputs to 
criterion outputs are computed. The bandwidth of the closed loop system will be approximately 
that of the open loop synthesis model. Figure 22 contains a block diagram of the synthesis model 
and the design weightings used for this example. 
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Select Regulated Variables: 
o Heading angle (Psi) 
o Sideslip angle (Beta) 

Define ideal Models: 

o Heading command processor: first order lag at 0.4 rad/sec 

- Heading command processor output is rate limited to limit 
bank angle at 30 degrees 

o Sideslip command processor: first order lag at 0.5 rad/sec 



Figure 21. Feedforward Controller Showing Regulated Variables and Ideal Models 


Psi command 



Synthesis Transmission Zeros: 

o Psi Criterion: 1 real zero: s = -3 

o Beta Criterion: 2 complex zeros: 
Wn = 2 Zeta = .8 


LQG Weighting Matrices: 



— — 


— — 

Q = 

2.5 

R = 

1 


2.5 


1 


Figure 22. LQG Synthesis Model Showing Criteria Outputs and Weighting Matrices 
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For this example, only the sideslip state was deemed insufficient for direct feedback. The approach 
taken was to use a complementary filter to estimate Beta. It was assumed that an air data sensor 
would give low frequency Beta while high frequency information could be derived from inertial 
data. A simple first order complementary filter at 0.1 rad/sec was proposed. Since the 
complementary filter does not affect the system stability it did not figure in the following analysis. 
A complete analysis of the performance of the Beta filter would require analysis of the sensor 
characteristics which is beyond this example. 

8.1.5 CLOSED LOOP HEADING CONTROLLER 

Figure 25 is a block diagram of the closed loop system for the heading control including the sideslip 
complementary filter. 


8.2 HEADING CONTROLLER RESULTS 

The following paragraphs present closed loop heading controller analysis results. The feedback 
controller gains were defined per the gain schedules presented in paragraph 8.1.4 using the structure 
given in Figure 25. Expanded analysis giving results for all 160 design flight conditions is presented 
in Plots 31 through 45 found in Appendix B. 

8.2.1 EIGENVALUES 

Figure 26 is a scatter plot of the eigenvalues of the closed loop heading controller. The dutch roll damping 
ratio, with a couple of minor exceptions, is above 0.6 for each of the fourteen design flight conditions. 

8.2.2 GAIN AND PHASE MARGINS 

Figure 27 shows the phase and gain margin characteristics for the rudder and aileron loops. Cross 
plots of the broken loop frequency response real and imaginary components show that no loops 
violate the region designating +/- 4dB and +/- 40 degrees. Similar results were achieved for the 
sensor loops. 

8.23 COVARIANCE RESPONSES 

Figure 28 shows the position covariance responses of aileron, rudder, heading, and sideslip to a unit 
magnitude lateral dryden gust. These data are plotted versus flight condition number. These 
responses, as well as the rate responses, meet the requirements. 
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Phi limiter 



Figure 25. Detailed Closed Loop Block Diagram of Heading Controller 


PSI-BETA CONTROLLER CLOSED-LOOP SCATTER PLOT 



Figure 26. Scatter Plot for Heading Controller Using Scheduled Feedback Gains 
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Figure 27. Aileron and Rudder Broken Loop Nyquist Plots for Heading Controller 



Figure 28. Covariance Responses of Heading Controller to 1 Ft / Sec Dryden Gust 


69 


( 


8.2.4 FREQUENCY RESPONSES 


Figure 29 shows the heading command to heading frequency responses for each of the fourteen flight 
conditions. Since the response is, by design, dominated by the ideal model, there is very little 
variation in response from condition to condition. 

8.2.5 TIME DOMAIN SIMULATION 

Figure 30 shows time history plots of aileron, rudder, heading, and sideslip in response to a step in 
heading command for condition 1. Plots 46-57 found in Appendix B present nonlinear simulation 
results for a 90 degree heading command change. Plots 46-51 show results for a low speed flight 
condition while Plots 52-57 show the same data for a high speed flight condition. 

8.3 LOCALIZER CONTROLLER DESIGN 

The localizer controller was designed using the heading controller developed in paragraph 8.1 as an 
inner loop. A classical root locus technique was used. The cross track error is sensed and fed back to 
command a change in heading. The inner loop heading controller then commands the airplane to fly 
to a new heading. 

8.3.1 CONTROLLER STRUCTURE AND GAINS 

The structure of the localizer controller is proportional plus integral cross track error fed back to form 
the heading command. Classical root locus was used to determine the localizer controller gains. 
(Figure 31 shows the block diagram of the localizer controller including the gain values.) The 
feedback gain is scaled by the inverse of airspeed to avoid aggressive control leading to overshoots 
during capture at high speed. 

83.2 LOCALIZER CAPTURE LOGIC 

When an autopilot approach for landing is made the airplane flies along under heading control 
until the localizer beam signal is received. Once the localizer signal is available, the plane must 
maneuver to capture the beam and follow it into the runway. The logic employed to transition from 
heading control to localizer track is described in the following. 
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PSI - BETA - A_ACT - R_ACT 


PSI CONTROLLER FREQUENCY RESPONSES FOR 14 CASES 
INPUT: PSI CMD 
OUTPUT: PSI 
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Figure 29. Frequency Responses of Heading Controller from psi Command to psi 


CLOSED-LOOP LATERAL STEP RESPONSE FOR CASE #001 
SCHEDULED GAINS 
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Figure 30. Time Response of Heading Controller for Heading Command Step 
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Figure 31. Block Diagram of Localizer Capture and Track Controller 



Figure 32. Outline of Localizer Capture Logic 
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Figure 32 outlines the localizer capture logic. The critical portion of the logic is testing to see if cross 
track error is less than a threshold defined by the velocity and maximum desired bank angle for autopilot 
turns. If the airplane is outside the threshold it turns toward the localizer beam until it is flying 
perpendicular to the beam with maximum closing speed. Once the airplane closes within the threshold 
it is commanded to turn parallel to the beam. When the difference between the airplane heading and the 
runway heading falls below 30 degrees, the localizer controller is engaged for final capture and track. 

8.4 LOCALIZER CONTROLLER RESULTS 

The following paragraphs present results from analysis of the closed loop localizer controller as 
defined by the block diagram in Figure 31. 

8.4.1 EIGENVALUES 

Figure 33 is a scatter plot of the eigenvalues of the closed loop localizer controller for each of the 
fourteen design flight conditions. Damping of 0.6 was achieved for all conditions. 

8.4.2 GAIN AND PHASE MARGINS 

Figure 34 shows the phase and gain margin characteristics for the rudder and aileron loops. Cross plots 
of the broken loop frequency response real and imaginary components show that no loops violate the 
region designating +/- 4dB and +/- 40 degrees. Similar results were achieved for the sensor loops. 

8.4.3 COVARIANCE RESPONSES 

Figure 35 sho ws the position covariance responses of aileron, rudder r heading 7 and sideslip to a unit 
magnitude lateral dryden gust. These data are plotted versus flight condition number. These 
responses, as well as the rate responses, meet the requirements. 

8.5 NONLINEAR SIMULATION OF LOCAUZER CAPTURE AND TRACK 

The time history plots in Figures 36 and 37 illustrate two localizer capture senarios. Figure 36 shows 
the ground track path followed for capture from an initial condition flying parallel to the localizer 
beam with a displacement of 15000 feet, while Figure 37 shows the ground track path followed for 
capture from an initial condition flying 90 degrees to the beam with initial displacement of zero. 

These two test cases illustrate the function of transition from heading control to localizer capture 
and track. 
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LOCALIZER CONTROLLER CLOSED-LOOP SCATTER PLOT 



Figure 33. Scatter Plot for Localizer Capture and Track Controller 
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Figure 34. Aileron and Rudder Broken Loop Nyquist Plots for Localizer Capture and Track Controller 
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LOCALIZER POSITION COVARIANCE RESPONSE TO 1 FT/SEC RMS 
DRYDEN TURBULENCE FOR CONDITION SPEED AND ALTITUDE 



Figure 35. Covariance Responses of Localizer Capture and Track Controller to 1 Ft / Sec Dryden Gust 


LOCALIZER CAPTURE FOR CASE #086 
30 DEGREES FLAP, LIGHTWEIGHT, 241 FT/SEC 
INITIAL CONDITION: Y - 15000. FT, PSI * 0. 
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Figure 36. Flight Path History for Capture from a Parallel Heading with 1500 Ft Initial Offset 
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LOCALIZER CAPTURE FOR CASE #086 
30 DEGREES FLAP. LIGHT WEIGHT, 241 FT/SEC 
INITIAL CONDITION: Y - 0. PSI - 90. 



X 
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Figure 37. Flight Path History for Capture from a Perpendicular Path with No Initial Offset 
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Additional non-linear simulation results for the localizer capture and track controller are presented 
in Plots 58-81 found in Appendix B. These 24 plots are divided into four groups (58-63, 70-75, and 76- 
81) corresponding to a low and high speed captures from a heading 45 degrees relative to the runway 
and low and high speed captures from a heading 90 degrees relative to the runway respectively. 
Each plot gives three traces: 1) still air, 2) a 20 ft/ sec crosswind flowing toward the runway 
centerline, and 3) a 20 ft/ sec crosswind blowing away from the runway centerline. A listing of the 
FORTRAN program used to generate these plots is found at the end of Appendix B. 

8.6 SUMMARY OF LATERAL AXIS DESIGN 

The integral LQG design process presented in this section is an efficient technique for the design of 
multiple input / multiple output control systems. Unlike classical root locus methods, requirements 
of more than one loop can be handled at a time. In addition, the systematic approach to selecting 
regulated variables and forming criteria outputs affords more insight than LQG techniques using 
state weighting only. The use of transmission zeros, both inherent in the plant and created by the 
designer, is crucial since they dictate the asymptotic nature of the closed loop eigenvalues. 

Although classical root locus techniques treating one loop at a time have been adequate in the past 
for the design of transport airplane control systems, the need is arising for multiple loop design 
methods. For example, in search of greater efficiency, airplane structural stiffness has been reduced 
leading to lower frequency flexible modes. The result is a challenging controls problem since control 
inputs now risk exciting structural modes. The technique presented here allows engineers with 
experience applying classical design methods to quickly learn a multiple loop design approach since 
the single loop compensation ideas they are familiar with are the same tools used to develop 
criterion outputs, the centraf feature of this method. 
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9.0 CONCLUSIONS 


An integrated autopilot/autothrottle control system has been developed. A two-degree-of-freedom 
approach was used to achieve a satisfactory design that offers the designer the fundamental 
advantages of: 

a . Simple mode switching logic 

b. Limiting done in feedforward loop, hence, no stability effects 

c. Uniform closed loop response throughout the flight envelope. 

Within the two-degree-of-freedom framework, the feedback regulator was designed using an 
integral LQR design technique, which offers a systematic approach to satisfy desired feedback 
performance requirements and guarantees stability margins in both control and sensor loops. 

The resulting feedback controller was discretized and implemented using a delta coordinate concept, 
which allows for transient free controller switching by initializing all controller states to zero and 
provides a simple solution for dealing with throttle limiting cases. 

In conclusion, it was shown, that a systematic top-down approach to complex control design problems 
combined with proper application of modern control synthesis techniques yields a satisfactory 
solution in a reasonably short time period. 
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APPENDIX A 


GLIDESLOPE MODEL COMPUTATIONS 

The computations shown here are based on [5]. 

From Figure 5.5: 

x = Hrad cot (GSA - GSE) 

He = x tan GSA 

He = H c - Hrad 

= x tan GSA - Hrad 

„ ( tan GSA \ 

RAD [tan (GSA - GSE) " 

and, using small angle approximations: 

„ u GSE 

H e Hrad gsa . GSE 


To compute He consider: 


He 


d TT dHc dx 
dt ** = dx ’ dt 


but 


dHc 

dx 


tan GSA = 


GSA 

57.3 


and 



H c = Vi 


GSA 

57.3 


FLARE PATH COMPUTATIONS 


Altitude command: 

Hc(x) = Ax^ + Bx^ i- C x + D 
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Altitude rate command: 


Vi , Vi Vi 

H c (x) = 3A — — x 2 + IB-*- + C — zr 
L xtd xtd X TD 

xjD = distance to touchdown from flare initiation = 1200 ft 

Vj = inertial speed of the aircraft 

Find A, B, C, D from boundary conditions: 

x = 0 x = 1200 

He = 45 ft He = 0 

H c = H H c = -2.5 fps 

D = 45 ft 


C = 


H/x = o xtd 

vi 


B = 2C - 3D - 45 


XTD 

Vi 


A = 


H xjD 
3Vi 


2B 

3 


C 

3 


SPEED HOLD COMPUTATIONS 


1 . Speed conversion 
MACH -> VTAS : 

By definition 

VTAS = MACH ■ c 


where 

c is a speed of sound 
VC AS -> VTAS I5j : 

VTAS = VC AS ^ — 

1 - 10' 5 H 

2. Limiter computations: 

V limiter : Consider E s (energy rate per pound of weight) equation: 


( 1 ) 


( 2 ) 
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(3) 


V 

E s = H + V • — 

S g 

Es _ H V T-D 
v " V + g _ w 

When thrust is constant: 


Es 

V 





H 
'g V 


H limit is 0.05 g = 1.6 fps^ 


Hence, 


VLIM 


32.2 -1.6 52 

V “ V 


(4) 


(5) 


( 6 ) 


V limiter : When throttles are at the limit, airplane's energy rate is constant, (4) , assuming 
drag doesn't change. Therefore, any speed change must be accomplished at the expense of 
climb rate. 

Let Hmin be the minimum climb rate desired. Then 


AE S _ Hmin - H VLIM = Q (7) 

V ‘ V + g 



In climb Hmin = 10 fps 


In descent Hmin = 0 fps 
In climb (max thrust): 



In descent (min thrust) 
H 

V min - g y 


( 8 ) 


(9) 


( 10 ) 
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NOMINAL FLIGHT CONDITION FOR CRUISE CONTROLLER DESIGN 


A MATRIX (5X5) 
-0. 422307E-01 
-0.279749E-02 
0 . 686186E-01 
0 .OOOOOOE+OO 
-0.196946B-06 
B MATRIX (5X5) 
-0. 661506E-01 
-0.321807E-01 
-0.297967E+01 
0 .OOOOOOE+OO 
O.OOOOOOE+OO 
C MATRIX (8X5) 
-0. 422307E-01 
0 . 100000E+01 
-0 . 196946E-0 6 
O.OOOOOOE+OO 
0.000000E+00 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
-0.122523E-02 
D MATRIX (8X5) 
-0.661506E-01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
-0.133329E-01 
STATES (5) 
'U' 'ALPHA' 
INPUTS (5) 
'DE' ' DSPL' 
OUTPUTS (8) 

' UD' 'UO' 


.STATE SPACE MODEL OF THE AIRPLANE 


-0.39I491E+00 
-0 . 542493E+00 
-0.207439E+01 
O.OOOOOOE+OO 
-0.132435E+02 

-0 . 507741E-02 
0. 942430E-02 
0.291736E-01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

-0 . 391491E+00 
O.OOOOOOE+OO 
-0 . 132435E+02 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
-0 . 223450E+00 

-0. 507741E-02 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0.386026E-02 

'Q' 'THETA' 


-0.775813E+00 
0 . 999048E+00 
-0 . 907272E+00 
0 . 100000E+01 
O.OOOOOOE+OO 

0.292118E-01 
-0.131920E-03 
0 .218543E-01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

-0.775813E+00 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 . 100000E+01 
O.OOOOOOE+OO 
-0 .179852E— 02 

0.292118E-01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 . 185995E-01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
-0.116222E-05 

'H' 


-0.560582E+00 

0.247901E-05 

0.170303E-03 

O.OOOOOOE+OO 

0.132435E+02 

0.422307E-01 
0.279749E-02 
-0.686186E-01 
O.OOOOOOE+OO 
0 . 19694 6E-06 

-0.560582E+00 
O.OOOOOOE+OO 
0 . 132435E+02 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 . 100000E+01 
-0 . 101611E— 02 

0.422307E-01 
—0 . 10000 0E+01 
0.196946E-06 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0.122523E-02 


'DTH' 'UG' ' ALFG' 

' HDO' 'HO' ' EPRO' 'QO' 'THO' 'NZ' 


END 

CONDITION# 

MACH 

FLAPS 

H 

GAMMA 

VTAS 

VCAS 

CG 

WEIGHT 

ALFA 


159 

0.7799900 

0.0000000 

35000.00 
-1.1285000E-05 

758.8500 

264.5300 

0.3000000 

110000.0 

3.349600 


Q 246.5400 

GEAR 0.0000000 

THRUST /THROTTLE 134.9190 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

0 . 391491E+00 
0 . 542493E+00 
0.207439E+01 
O.OOOOOOE+OO 
0 . 132435E+02 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

0.100000E+01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

0 . 391491E+00 
O.OOOOOOE+OO 
0 . 132435E+02 
JO. OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0.223450E+00 


NOMINAL FLIGHT CONDITION FOR LANDING CONTROLLER DESIGN 


A MATRIX (5X5) 
-0.4 60346E-01 
-0 . 736570E-01 
-0 . 100121E-02 
O.OOOOOOE+OO 
-0.3 67 103E-0 7 
B MATRIX (5X5) 
0 . 1 942 68E-02 
-0 . 519616E-01 
-0 . 104851E+01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
C MATRIX (7X5) 
-0 . 460346E-01 
0.100000E+01 
-G.367103E-07 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
D MATRIX (7X5) 
0 . 1942 68E-02 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
STATES (5) 
'U' 'ALPHA' 
INPUTS (5) 
'DE' 'DSPL' 
OUTPUTS (7) 
'UD' 'UO' 


.STATE SPACE 
0 . 383599E+00 
—0.7 55050E+00 
-0 . 339359E+00 
O.OOOOOOE+OO 
-0.394695E+01 

-0.971315E-02 
0 . 343995E-01 
0 . 344566E-01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

0 . 383599E+00 
O.OOOOOOE+OO 
-0.394695E+01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

-0 . 971315E-02 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

'Q' 'THETA' 


MODEL OF THE 1 
—0 . 374943E— 01 
0 . 100284E + 01 
-0.739641E+00 
0 . 100000E+01 
O.OOOOOOE+OO 

0.224840E+00 
-0.554790E-03 
0 . 12 8302E+00 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

-0.374943E-01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 . 100000E+01 
O.OOOOOOE+OO 

0.224840E+00 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

0.226125E-01 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

'H' 


-0 . 5614 53E+00 
-0.722729E— 04 
-0 .3010 61E-03 
O.OOOOOOE+OO 
0.394695E+01 

0 . 46034 6E-01 
0.736570E-01 
0.100121E-02 
O.OOOOOOE+OO 
0 . 367103E-07 

-0.561453E+00 
O.OOOOOOE+OO 
0.394695E+01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 . 100000E+01 

0 . 460346E-01 
-0.100000E+01 
0 . 367103E-07 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 


'DTH' 'UG' 'ALFG' 

'HDO' 'HO' 'EPRO' 'QO' 'THO' 


END 

CONDITION# 

MACH 

FLAPS 

H 

GAMMA 

VTAS 

VCAS 

CG 

WEIGHT 

ALFA 

Q 

GEAR 

THRUST /THROTTLE 


102 

0.2025700 
39.99900 
100.0000 
-2 . 1035000E-06 
226.1600 

134.0000 
0.3000000 

80000.00 
0.5426800 
61.41800 
1.000000 
559.0780 


O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

O.OOOOOOE+OO 

-0 . 383599E+00 
0 . 755050E+00 
0 . 339359E+00 
O.OOOOOOE+OO 
0.394695E+01 

O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 . 100000E+01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 

-0 . 383599E+00 
O.OOOOOOE+OO 
0 . 394695E+01 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 


ELEVATOR GAINS FOR LANDING CONTROLLER 



A + B ( Q c- Q c o ) 

1 + C(Q c -Q Co ) 


K = 1.2 


(0.841-0. 0023(Q C - 43. 5)) 
K u= 1+0.00223 (Q c - 43.5) 

_ D + E (g-a 0 ) 

H 1 + F (a - a Q ) 

3. 44 + 0. 0229(Q C - 43. 5) 
K Q = 1 + 0. 0183(Q c -43.5) 

K e = 7.1 

g+_h(q c -Qc o) 

EPR 1 + I /Q - Q_ \ 


H + K(a-a 0 ) 
K h “ 1 + L (a - a Q ) 

Q_ = NOMINAL DYNAMIC PRESSURE 
c o 


WHERE A, B, C, D, E, F, G, H, I, ], K, L a function of flap position 
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THROTTLE GAINS FOR LANDING CONTROLLER 


A + B ( Q c " Q c 0 ) 

1 +C ( Q c - Q c 0 ) 

D + E ( Q c - Q C J 

1 +F (° C " °c 0 ) 

- 2. 1 53 + 0.0078(Q C - 43.5) 
K u = 1.0 - 0.0033(Q C - 4 3.5) 


G * H (Q c - Q Cp ) 
H = 1 + I(Q C - Q Cq ) 


I + K(Q C - Q Cq) 
Q = 1 + L (°‘ - °c 0 ) 


0.422 - 0. 238(a + 2.81) 
9 ~ 1.0 + 0.007(a + 2.81) 


- 183 + 0. 22CTHTRM - 13.87) 
EPR 1.0 - 0. 0046(THTRM - 13.87) 


-yc 

« - 1 + P(Q C - Q Co ) 


K 


K 


H 


WHERE A, B, C, D, E, F, G, H, I, J, K, L, M, N, P are a function of flap position 
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DISCRETIZED CONTROLLER COMPUTATIONS 


1. Feedback controller: 


A 


d 


-A -A. 

e * 0 0 

1 0 

0 0 1 



B 


d 




0 


0 

AT 

0 


0 

0 

1 


2. Feedforward controller: 


A 


d 


.9851 - 1.117E-3 

4.96E-2 1 


B 


d 


1. 117E- 3 
28E-5 
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LINEAR ANALYSIS RESULTS 
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Plot 1. Open Loop Short Period Poles Scatter Plot for Conditions 113 - 160 



Plot 2. Open Loop Short Phugoid Poles Scatter Plot for Conditions 113 - 160 


IMAG 



Plot 4. Cruise Controller Closed Loop Phugoid Eigenvalues Scatter Plot 
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COND 

Plot 5. Cruise Controller Broken Loop Analysis Plot of Gain Margins 
vs Flight Condition for All Control and Sensor Loops 
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Plot 6. Cruise Controller Broken Loop Analysis Plot of Phase Margins 
vs Flight Condition for All Control and Sensor Loops 
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Plot 14- Cruise Control Law Turbulence Response Analysis Covariance Response 
to Dryden Turbulence vs Flight Condition (Altitude in Ft) 
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Plot 19. Cruise Control Law Turbi 
to Dryden Turbulence vs 
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NONLINEAR TESTING RESULTS 


ALTITUDE FROM TRIM, FT 



AIRSPEED FROM TRIM, FT/SEC 


Z 0.006 



Plot 22. Nonlinear Analysis: Glideslope Capture Airplane is Trimmed at 2000 Ft, Flaps 25 



Plot 23. Nonlinear Analysis: Glideslope Capture Airplane is Trimmed at 2000 Ft, Flaps 25 
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Plot 24. Nonlinear Analysis: Glideslope Capture Airplane is Trimmed at 2000 Ft, Flaps 25 



TIME 


SOE: SPEED ON ELEVATOR 

Plot 25. Nonlinear Analysis: Cruise Controller Response to Simultaneous Altitude and 

Speed Commands with Throttle Limiting (Nominal Plant is Trimmed at 35,000 Ft) 




TIME 

Plot 26. Typical Submode Logic Events for Flight Level Change Mode of Conventional Autopilot / Autothrottle 


2 

S 


o 

DC 

U. 

c/> 

z 

o 

H 

< 

E 

$ 


50 

45. 

40. 

35. 

30. 

25. 

20 . 

15. 

10 . 

5. 

0 . 

-5._ 

-50 


THROTTLE POSITION IN DEG 

♦ 



V.STEP 


50 150 250 350 

TIME 


Plot 27. Nonlinear Analysis: Cruise Controller Response to Slimultaneous Altitude and 
Speed Commands with Throttle Limiting (Nominal Plant is Trimmed at 35,000 Ft) 
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APPENDIX B 


LATERAL AXIS RESULTS: 

AILERON AND RUDDER GAIN SCHEDULES FOR CRUISE CONTROLLER 
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Plot 30A. Design Aileron Feedback Gains for Each Flight Condition 
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Plot 30B. Design Rudder Feedback Gains for Each Flight Condition 
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Rudder Gain Schedules For Cruise Controller 


















LINEAR ANALYSIS RESULTS 
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Plot 31. Open Loop Poles for Flight Conditions 1-80 
FOUR STATES: P, R f PHI, BETA 
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Plot 32. Open Loop Poles for Flight Conditions 81 - 160 
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Plot 33. Closed Loop Poles for Flight Conditions 1-80 



Plot 34. Closed Loop Poles for Flight Conditions 81 - 160 
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Plot 35. Broken Loop Phase Margins for Aileron and Rudder Loops at All Flight Conditions 
150. 


100 


50 


-50 


-100 


-150 


-200 


.sa$ 

ft? 


?P° d D 

] n n nn 

tct-^ 

3 °°n° 

5gj 

S 

T3 

2^, 





















! 

S&i 




D uu 




3qQ &S 
□ □ 
D° 

O d dDQ 
□ a .J 

]DOd 

5> d od” 

n aD ^Ql 

j 






d w u 
3 

□ 

□ 

□ 

CT 


□ 

* 


40 


80 

COND 


120 160 

Plot 36. Broken Loop Phase Margins for All Sensor Loops at All Flight Conditions 
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Plot 37. Broken Loop Gain Margins for All Sensor Loops at All Flight Conditions 
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Plot 38. Aileron and Rudder Broken Loop Bandwidths at All Flight Conditions 
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Plot 39. Sensor Broken Loop Bandwidths at All Flight Conditions 
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Plot 40. 


Closed Loop Command Frequency Response From ycMD To y 
(Typical of all conditions) 
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Plot 41. Closed Loop Adverse Command Frequency Response From \|/CMD To (3 
(Typical of all conditions) 


SCH COMP(.I) ACT(15,15) INPUT: VGUST 



Plot 42. Position Covariance Responses of V(/CMD To P to 1 ft/sec RMS 
Dryden Turbulence at Each Flight Condition 
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5. Rate Covariance R esponse of Aileron and Rudder to 1 ft/sec RMS 
Dryden Turbulence at Each Flight Condition 
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TIME SIMULATION RESULTS FOR 


HEADING CONTROLLER 
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TIME SIMULATION RESULTS FOR 
LOCALIZER CAPTURE AND TRACK CONTROLLER 

.1 4 
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„ „ LOCALIZER CAPTURE FOR CASE #001 

THREE CASES: STILL AIR, 20 FT/SEC CROSSWIND, -20 FT/SEC CROSSWIND 
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Plot 58. Ground Track Response During Localizer Capture and Track 

Initial condition: heading 45° from runway heading 
Low speed flight condition: flaps = 1, VCAS « 138 knots 

Tuner ^ __L° CA Jt ,zer CAPTURE FOR CASE #001 
THRE /t££§§§l STILL AIR 20 FT/SEC CROSSWIND, -20 FT/SEC CROSSWIND 
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Plot 59. Ground Track Heading Response During Localizer Capture and Track 

Initial condition: heading 45° from runway heading 
Low speed flight condition: flaps = 1, VCAS = 138 knots 
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PHI 




Initial condition: heading 45° from runway heading 
Low speed flight condition: flaps = 1, VCAS = 138 knots 
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Initial condition: heading 45° from runway heading 
Low speed flight condition: flaps = 1, VCAS = 138 knots 
LOCALIZER CAPTURE FOR CASE #001 


THREE CASES: STILL AIR. 20 FT/SEC CROSSWIND. -20 FT/SEC CROSSWIND 
(POSITIVE CROSSWI^D BLOWS AIRPLANE AWAY FROM RUNWAY) 



Plot 63. Rudder Response During Localizer Capture and Track 

Initial condition: heading 45° from runway heading 
Low speed flight condition: flaps = 1, VCAS =138 knots 
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Initial condition: heading 45° from runway heading 
High speed flight condition: flaps = 1, Vcas = 210 knots 
LOCALIZER CAPTURE FOR CASE #004 

THREE CASES: STILL AIR. 20 FT/SEC CROSSWIND, -20 FT/SEC CROSSWIND 



Plot 65. Ground Track Heading Response During Localizer Capture and Track 

Initial condition: heading 45° from runway heading 
High speed flight condition: flaps = 1 , Vcas = 210 knots 
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High speed flight condition: flaps = 1, Vcas = 210 knots 



Initial condition: heading 45° from runway heading 
High speed flight condition: flaps = 1, Vcas = 210 knots 
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LOCALIZER CAPTURE FOR CASE #004 
THREE CASES: STILL AIR, 20 FT/SEC CROSSWIND, -20 FT/SEC CROSSWIND 
(POSITIVE CROSSWIND BLOWS AIRPLANE AWAY FROM RUNWAY} 



Initial condition: heading 45° from runway heading 
High speed flight condition: flaps * 1, Vcas = 210 knots 
LOCALIZER CAPTURE FOR CASE #004 .... 


THREE CASES: STILL AIR. 20 FT/SEC CROSSWIND, , -20 FT/SEC CROSSWIND 
(POSITIVE CROSSWlfoD BLOWS AIRPLANE AWAY FROM RUNWAY) 



Plot 69. Rudder Response During Localizer Capture and Track 

Initial condition: heading 45° from runway heading 
High speed flight condition: flaps = 1, Vcas = 210 knots 
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LOCALIZER CAPTURE FOR CASE #001 

THREE CASES: STILL AIR, 20 FT/SEC CROSSWIND, -20 FT/SEC CROSSWIND 
(P^TIVE CROSSWIND BL6WS AIRPLANE AWAY FROM RUNWAY CENTERLINE) 



Initial condition: 15000 ft offset from runway 

Low speed flight condition: flaps = 1, Vcas = 138 knots 


LOCALIZER CAPTURE FOR CASE #001 

THREE CASES: STILL AIR. 20 FT/SEC CROSSWIND, -20 FT/SEC CROSSWIND 
(PO^TIVE CROSSWIND BL6WS AIRPLANE AWAY FROM RUNWAY CENTERLINE) 



Plot 71. Ground Track Heading Response During Localizer Capture and Track 
Initial condition: 15000 ft offset from runway 
Low speed flight condition: flaps = 1, Vcas = 138 knots 
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LOCALIZER CAPTURE FOR CASE #001 



Initial condition: 15000 ft offset from runway 
Low speed flight condition: flaps = 1, Vcas = 138 knots 



^ * o jt 

Initial condition: 15000 ft offset from runway 
Low speed flight condition: flaps = 1, Vcas * 138 knots 
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Initial condition: 15000 ft offset from runway 
Low speed flight condition: flaps = 1, Vcas = 138 knots 



Plot 75. Rudder Response During Localizer Capture and Track 
Initial condition: 15000 ft offset from runway 
Low speed flight condition: flaps = 1, Vcas = 138 knots 
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Initial condition: 15000 ft offset from runway 

High speed flight condition; flaps = 1, Vcas = 210 knots 


LOCALIZER CAPTURE FOR CASE #004 

THREE CASES: STILL AIR, 20 FT/SEC CROSSWIND. -20 FT/SEC CROSSWIND 
(POSITIVE CROSSWl^D BLOWS AIRPLANE AWAY FROM RUNWAY) 



Plot 77. Ground Track Heading Response During Localizer Capture and Track 
Initial condition: 15000 ft offset from runway 
High speed flight condition: flaps = 1, Vcas = 210 knots 
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LOCALIZER CAPTURE FOR CASE #004 

THREE CASES: STILL AIR, 20 FT/SEC CROSSWIND, -20 FT/SEC CROSSWIND 
(POSITIVE CROSSWIND BLOWS AIRPLANE AWAY FROM RUNWAY) 



Initial condition: 15000 ft offset from runway 

High speed flight condition: flaps * 1, Vcas = 210 knots 


LOCALIZER CAPTURE FOR CASE #004 

THREE CASES: STILL AIR. 20 FT/SEC CROSSWIND. -20 FT/SEC CROSSWIND 
(POSITIVE CROSSWlKlD BLOWS AIRPLANE AWAY FROM RUNWAY} 



Plot 79. Bank Angle Response During Localizer Capture and Track 
Initial condition: 15000 ft offset from runway 
High speed flight condition: flaps * 1, Vcas = 210 knots 
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LOCALIZER CAPTURE FOR CASE #004 



Plot 80. Aileron Response During Localizer Capture and Track 
Initial condition: 15000 ft offset from runway 
High speed flight condition: flaps * 1, Vcas = 210 knots 



Initial condition: 15000 ft offset from runway 

High speed flight condition: flaps = 1, Vcas * 210 knots 
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SIMULATION PROGRAM 
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PROGRAM SIM 


C 

C This program performs a simulation of the lateral axis of a 737-200. 

C Two modes of operation are available: 

C 1) Heading / Sideslip control 

C 2) Localizer capture and track / Sideslip control 

C Note that mode 1 is used as the inner loop for mode 2. 

C 

C A linear point simulation plant model is used along with a linear controller 
C whose gains are scheduled as functions of the trim flight condition. 

C 

C 

CHARACTER LINE*80 

REAL A(4,4), B (4 , 3) , MACH, G(2,6), X(6), U_C(2>, U__X<3), 

+ XA_DOT ( 4 ) , XB_DOT ( 4 ) 

C 

C I/O files: 

C Unit 5: Linear plant model and trim flight condition data 

C Unit 6: Time simulation output data 

C 

OPEN (UNIT-5, FILE-' OL/LATOO 4 .MDL' , STATUS-' OLD' ) 

OPEN (UNIT— 6, FILE-' SIM. GGP' , STATUS-' UNKNOWN' ) 

C 

C Read in plant model and flight condition data for gain scheduling 
C 

100 READ (5/ 8010) LINE 

IF (LINE (1:1) .NE. 'A') GO TO 100 
DO 110 I - 1,4 

READ (5, *) (A (I, J) , J-1,4) 

110 CONTINUE 
C 

200 READ (5,8010) LINE 

IF (LINE (1:1) .NE. 'B') GO TO 200 
DO 210 1-1,4 

READ (5, *) (B(I,J), J— 1, 2 ) 

210 CONTINUE 
C 

300 READ (5,8010) LINE 

IF (LINE (6: 8) .NE. 'DFM' ) GO TO 300 
READ (LINE(14:23) ,8020) FLAP 
C 

400 READ (5,8010) LINE 

IF (LINE (6: 9) .NE. 'MACH') GO TO 400 
READ (LINE(14:23) , 8020) MACH 
C 

500 READ (5,8010) LINE 

IF (LINE (6:7) .NE. 'VP') GO TO 500 
READ (LINE (14 :23) , 8020) VTAS 
C 

600 READ (5,8010) LINE 

IF (LINE (6: 9) .NE . 'VCAS') GO TO 600 
READ (LINE (14:23) ,8020) VCAS 
C 
C 

C Build lateral gust input 
C 

DO 700 1-1,4 

B(I,3) - - 57.3 * A(I, 1) / VTAS 

700 CONTINUE 
C 
C 

C Aileron feedback gains: scheduled as functions of trim flight condition 


c 

G (1, 1) 

- Aileron 

sideslip gain 

c 

G (1, 2) 

- Aileron 

roll rate gain 

c 

G (1, 3) 

= Aileron 

bank angle gain 

c 

G ( 1 , 4) 

= Aileron 

yaw rate gain 

c 

G ( 1 , 5) 

= Aileron 

integral sideslip error gain 

c 

G(l, 6) 

- Aileron 

heading error gain 


C 

G (1, 1) = -4.2 


C 

IF (VCAS -LT. 125.) THEN 
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G ( 1 , 2 ) - -.80 

ELSE IF (VCAS .LT. 200.) THEN 

G ( 1, 2) - -.80 + {VCAS - 125.) • .42 / 75. 

ELSE IF (VCAS . LT . 360.) THEN 

G(l,2) - -.38 + (VCAS - 200.) * .20 / 160. 

ELSE 

G (1, 2) - -.18 

END IF 
C 

IF (MACH .LT. 0.15) THEN 
G (1, 3) - -1.4 

ELSE IF (MACH . LT . 0.31) THEN 

G (1, 3) - -1.4 + (MACH - 0.15) * 0.8 / 0.16 

ELSE IF (MACH .LT. 0.9) THEN 

G (1, 3) * -0.6 + (MACH - 0.31) * 0.3 / 0.59 

ELSE 

G (1, 3) = -0.3 

ENDIF 
C 

IF (MACH .LT. 0.1) THEN 
G<1,4) - -1.9 

ELSE IF (MACH .LT. 0.36) THEN 

G (1, 4) - -1.9 + (MACH - 0.1) * 1.25 / 0.26 

ELSE IF (MACH .LT. 0.9) THEN 

G (1, 4) - -0.65 + (MACH - 0.36) * 0.275 / 0.54 

ELSE 

G (1, 4) - -0.375 

ENDIF 
C 

G (1/ 5) » 0.0 

C 

G (1, 6) - -4.1 

C 

c 

C Rudder feedback gains: scheduled as functions of trim flight condition 

C G (2 , 1) - Rudder sideslip gain 

C G (2, 2) - Rudder roll rate gain 

C G (2, 3) * Rudder bank angle gain 

C G (2 , 4 ) - Rudder yaw rate gain 

C G (2, 5) * Rudder integral sideslip error gain 

C G(2,6) - Rudder heading error gain 

C 

G (2 , 1) * -4.1 - FLAP * 1.5 / 40. 

C 

G (2, 2 ) - .1625 * VCAS / 100. - .5625 

IF (G(2,2) .LT. -0.4) G(2,2) * -0.4 

IF (G (2/ 2) .GT. -0.075) G(2,2) - -0.075 

C 

G(2, 3) - G <2 , 2) 

C 

IF (VCAS .LT. 100.) THEN 
G<2,4) = 5.0 

ELSE IF (VCAS .LT. 200.) THEN 

G (2 f 4) - 5.0 - (VCAS - 100.) * 1.6 / 100 . 

ELSE IF (VCAS .LT. 360.) THEN 

G (2/ 4) - 3.4 - (VCAS - 200.) * .65 / 160. 

ELSE 

G (2, 4) - 2.75 

ENDIF 
C 

G (2, 5) - -5.5 

C 

G (2 , 6) - 2.2 


Initail conditions for simulation 

Actuators 

RUD_C * 0 . 

RUD_X * 0 . 

AIL_C = 0 . 
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c 

c 

c 


AIL__X 

States 


BETA 

* 

0. 

P 


0. 

PHI 

* 

0. 

R 

* 

0. 

PSI 

«= 

0. 

Y POS 

SB 

15000 

X POS 


0. 


C 

c 

c 


c 

c 

c 


Beta filter and integrator 


BETA_F 

BETA_I 

Commands 

PSI_C 

PSI_X 

BETA_C 

BETA_X 

Limits 

PHI_MAX 
R_MAX 
P__MAX 
RDJMAX 
Y THRSH 


0 . 

0 . 


30. 
32.2 
7 . 

32.2 

57.3 


* PHI_MAX / VTAS 

* P_MAX / VTAS 

* (VTAS +75.) * (VTAS +75.) / 

(32.2 * PHI MAX) 


Simulation controls 


I LOCAL 
T_INIT 
T_FINAL 
DT 

NPRINT 

IPRINT 


1 

- 1 . 

150. 

0.05 

20 

0 


C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


WRITE ( 6 , 9100) 

Simulation loop 

DO 5000 T * T_INIT, T__F INAL , DT 
Set disturbance levels 


X_WIND = 0 . 

Y__WIND * T 

IF (YJWIND .LT. 0.) 
IF (Y WIND .GT. 20.) 


Y_WIND 
Y WIND 


0 . 
20 . 


Translate earth fixed wind into body coordinate lateral gust 

V_GUST * - X_WIND * SIN (PS 1/57 . 3 ) + Y_WIND * COS (PSI/57 . 3 ) 

Set heading command: 

I LOCAL = 0 gives heading controller 
ILOCAL = 1 gives localizer controller 

IF (ILOCAL .EQ. 0) THEN 

IF (T .GE. 0 . ) PSI_C - 90 . 

ELSE 
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no o o non non non non o oooooooooooo 


YABS - ABS ( Y__POS ) 

YSIGN - 1. 

IF (YABS .GT. 1.) YSIGN - YABS / Y_POS 
PSI__ABS “ ABS (PSI) 

Check to see if cross track is greater than threshold value 
If greater, command turn toward runway 
If less, 

a) command heading 30 degrees from runway if heading toward runway 

b) command heading 5 degrees from runway if heading away from runway 

c) monitor for switching to localizer controller when: 

i) localizer controller turn rate greater than 

heading controller turn rate 

ii) heading is less than 35 degrees 

iii) heading is in direction of the runway centerline 

IF (YABS .GT. YJTHRSH) THEN 

PSI_C - PSI_C + 10. * (-YSIGN*90. - PSI_C) * DT 

LATCH » 0 
ELSE 

P S I_DOTA - 10. * (— YSIGN*30 . - PSI_C) 

IF (PS I* YSIGN. GT . 0 . ) PSI_DOTA - 10. * <-YSIGN*5. - PSI_C) 

C_DOT - 0.04 * Y_POS + VTAS * SIN (PSIGT/57 . 3) 

PSI_DOTL - (-0.1 * 57.3 / VTAS) * C_DOT 

IF ( ( (ABS (PSI_DOTL) .GT. ABS <PSI_DOTA) ) .AND. 

+ (ABS(PSI_X) .LT. 35.) .AND. (PSI_X*YSIGN .LT. 0.)) 

+ .OR. (LATCH .EQ. 1) ) THEN 

PSI_C - PSI_C + PSI DOTL * DT 

LATCH - 1 
ELSE 

PSI_C - PSIJ3 + PSI_DOTA * DT 

ENDIF 
ENDIF 

ENDIF 


Apply roll rate and roll angle limits to heading command 


PSI_R_MAX - 
PSI_R_MIN - 
PSI_RATE = 

IF (PSI_RATE 
IF (PSI_RATE 
IF (PSI_RATE 
IF (PSI_RATE 
PSI X = PSI X 


PS I — RATE + RD_MAX * DT 
PSI__RATE - RD_MAX • DT 
0.4 * (PSI_C - PSI_X) 

. GT . PS I_R_MAX ) P S I_RATE 
. LT . PSI_R_MIN) PSI_RATE 
. GT . R_MAX) P S I_RATE 

. LT . -R__MAX) PS I_RATE 

PSI RATE * DT 


PSI_R_MAX 
PSI_R_MIN 
R__MAX 
-R MAX 


Compute sideslip command filter output 


BETA_X - BETA__X + 0.5 * (BETA_C - BETA_X) * DT 


Form heading and sideslip errors and integrate sideslip error 


PS I__ERR - PSIGT - PSI_X 

BETA__ERR = BETA_F - BETA_X 
BETA I * BETA I + BETA ERR * DT 


Compute aileron and rudder commands using feedback signals and gains 

X(l) * BETA_F 
X(2) * P 

X (3 ) « PHI 

X ( 4 ) - R 

X ( 5) - BETA_I 

X ( 6) = PSI_ERR 

CALL PRODUCT (G, 2, 6, X, VJZ) 

AIL_X * AIL_X + 15. * (U__C(1) - AILJC) * DT 
RUD_JX = RUD_X + 15. * (U_C(2) - RUD X) * DT 


Compute plant model response to control inputs and wind disturbance 
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xu) - 

BETA 





U X(l) 

- 

AIL X 





U X (2 ) 

« 

RUD X 





U X (3) 

* 

V GUST 





CALL PRODUCT (A, 4, 4 

, X, XA DOT) 




CALL PRODUCT (B, 4, 3 

, U_X, XB_DOT) 




BETA 


BETA + 

(XA DOT ( 1) + XB 

DOT ( 1 ) ) 

★ 

DT 

P 

wm 

P + 

(XA DOT (2) + XB* 

“DOT (2) ) 

* 

DT 

PHI 

* 

PHI + 

(XA DOT (3) + XB" 

DOT (3) ) 

★ 

DT 

R 

* 

R + 

(XA DOT (4) + XB‘ 

"DOT (4) ) 

★ 

DT 

PSI 

- 

PSI + 

R * DT 




B AIR 

* 

BETA - 

(57.3 / VTAS) * ^ 

7_GUST 



PSI AIR 

s 

PSI + B 

AIR 




X VEL 

= 

VTAS * COS (PSI AIR/57.3) 

+ X WIND 


Y VEL 

ar 

VTAS * SIN (PSI AIR/57.3) 

+ Y WIND 


X POS 

- 

X POS + 

X VEL * DT 




Y POS 


Y POS + 

Y VEL * DT 




PSIGT 

at 

57.3 * ATAN2 (Y_VEL,X_VEL) 





c 

c 

c 

c 


c 

c 

c 


Complementary filter for sideslip estimate 
air data sideslip - break frequency at 0, 


using initerial 
1 radians 


sideslip rate and 


B _ D0T - 1 
BETA F 


XAJDOT(l) 
BETA F + 


+ XB_DOT(l) 
(B DOT I + 0 , 


1* (B AIR - BETA F) ) * DT 


Write data to output file for post processing 


IF ((IPRINT .EQ. 0) .OR. (IPRINT 
WRITE (6,9200) T, BETA, 

+ PSI, PSIGT, 

+ PSI_C, PSI_X, 

+ B_AIR, BJDOT_I, 

+ AIL_X, RUD_X, 

IPRINT - 1 
ELSE 

IPRINT * IPRINT + 1 
ENDIF 


.GE. 


NPRINT) ) 
P, 


X_POS, 

BETA__C, 

VJ3UST, 

BETA_I, 


THEN 
PHI, 
Y_POS, 
BETA_X, 
X_WIND, 
I LOCAL, 


R, 

YJTHRSH, 

BETA_F, 

Y_WIND, 

LATCH 


C 

c 

5000 

C 


CONTINUE 


End of simulation loop 
WRITE (6, 9300) 


C 

C 

8010 

C 

8020 

C 

C 

9100 

C 


FORMAT (A80) 
FORMAT (G1 0.5) 


FORMAT ('$ LOCALIZER CAPTURE FOR CASE #004',/, 
+ '$ NO CROSSWIND (STILL AIR)',/, 

+ '$ ',/, 

+ '$ 20 FT/SEC CROSSWIND 
+ '$ (POSITIVE CROSSWIND 
+ 'SIM',/, 


(RAMPED IN DURING FIRST 20 SECONDS)',/, 
BLOWS AIRPLANE AWAY FROM RUNWAY)',/, 




C 

9200 


+ 

r 

T 

BETA 

P 

PHI 

R 

MORE 

+ 

9 

PSI 

PSIGT 

X POS 

Y POS 

Y THRSH 

MORE 

+ 

f 

PSI c 

PSI X 

BETA C 

BETA X 

BETA F 

MORE 

+ 

f 

B AIR 

B DOT I 

V GUST 

X WIND 

Y WIND 

MORE 

+ 

9 

AIL X 

RUD X 

BETA I 

I LOCAL 

LATCH' ) 



FORMAT ( IX, 5G15 . 5 , / , 

+ 2X, 5G15 . 5, / , 

+ 2X, 5G15 . 5, /, 

+ 2X, 5G15 . 5, /, 

+ 2X, 3G15 . 5, 2115) 
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9300 FORMAT (' *EOF' ) 

C 

END 

C 

c 

C Subroutine to form product of matrix and vector 
C 

SUBROUTINE PRODUCT (A, NROW, NCOL, X, Y) 

C 

REAL A (NROW, NCOL), X(NCOL), Y (NROW) 

C 

C 

DO 200 I * 1, NROW 
Y(I) - 0. 

DO 100 J - 1 , NCOL 

Y(I) - Y ( I) + A (I, J) * X ( J) 

100 CONTINUE 

200 CONTINUE 
C 
C 

END 
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